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We consider n x n real symmetric and Hermitian Wigner ran- 
dom matrices n~^^^W with independent (modulo symmetry condi- 
tion) entries and the (null) sample covariance matrices n~^X*X with 
independent entries of m x n matrix X. Assuming first that the 4th 
cumulant (excess) K4 of entries of W and X is zero and that their 
4th moments satisfy a Lindeberg type condition, we prove that linear 
statistics of eigenvalues of the above matrices satisfy the central limit 
theorem (CLT) as n ^ cx), m 00, m/n ^ c G [0, 00) with the same 
variance as for Gaussian matrices if the test functions of statistics are 
smooth enough (essentially of the class C^). This is done by using 
a simple "interpolation trick" from the known results for the Gaus- 
sian matrices and the integration by parts, presented in the form of 
certain differentiation formulas. Then, by using a more elaborated 
version of the techniques, we prove the CLT in the case of nonzero 
excess of entries again for essentially test function. Here the vari- 
ance of statistics contains an additional term proportional to K4. The 
proofs of all limit theorems follow essentially the same scheme. 



1. Introduction. The central limit theorem (CLT) is an important and 
widely used ingredient of asymptotic description of stochastic objects. In the 
random matrix theory, more precisely, in its part that deals with asymptotic 
distribution of eigenvalues {a["^}JL]^ of random matrices of large size n, 
natural objects to study are linear eigenvalue statistics, defined via test 
functions : R — > C as 
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A. LYTOVA AND L. PASTUR 



The question of fluctuations of linear eigenvalue statistics of random matrices 
was first addressed by Arharov [3] , who announced the convergence in prob- 
ability of any finite collection of properly normalized traces of powers of sam- 
ple covariance matrices in the case where the numbers of rows and columns 
of the data matrix are of the same order [see formulas (4.1), (4.2) and (4.7) 
below]. The result was restated and proved by Jonsson [16]. However, the 
explicit form of the variance of the limiting Gaussian law was not given in [3] 
and [16]. In 1975 Girko considered the CLT for the traces of resolvent of the 
Wigner and the sample covariance matrices by combining the Stieltjes trans- 
form and the martingale techniques (see [12] for results and references). In 
particular, an expression for the variance of the limiting Gaussian laws was 
given, although the expression is much less explicit than our formulas (3.92) 
and (4.65) for ip{X) = (A — z)~^ , ^ 0. In the last decade a number of re- 
sults on the CLT for linear eigenvalue statistics of various classes of random 
matrices has been obtained (see [2, 5, 7, 8, 10, 11, 13, 15, 16, 17, 18, 28, 29, 30] 
and [31] and the references therein). 

A rather unusual property of linear eigenvalue statistics is that their vari- 
ance remains bounded as n — > oo for test functions with bounded derivative. 
This has to be compared with the case of linear statistics of independent 
and identically distributed random variables l^;"^}]^^, where the variance is 
linear in n for any bounded test function. This fact is an important element 
of the ideas and techniques of the proof of the CLT for 

(1.2) {MnM - E{N^y]})/{Vs.r{KM})'^\ 

viewed as a result of addition of large numbers of small terms (see, e.g., 
[14], Chapter 18). On the other hand, since the variance of linear statistics 
of eigenvalues of many random matrices is bounded in n, the CLT, if any, 
has to be valid for statistics (1.1) themselves (i.e., without any n-dependent 
normalizing factor in front), resulting from a rather subtle cancelation be- 
tween the terms of the sum. One can also imagine that the cancelation is not 
always the case, and indeed it was shown in [24] that the CLT is not neces- 
sarily valid for so-called Hermitian matrix models, for which non-Gaussian 
limiting laws appear in certain cases even for real analytic test functions. 

In this paper we prove the CLT for linear eigenvalue statistics of two 
classes of random matrices: the Wigner matrices n~^^'^W, where W are nxn 
real symmetric random matrices with independent (modulo symmetry con- 
ditions) entries (typically n-independent) and the matrices n~^X^X, where 
X are m x n matrices with independent (and also typically n-independent) 
entries. We will refer to these matrices as the Wigner and the sample covari- 
ance matrices, respectively. The case, where the entries of W are Gaussian 
and the probability law of W is orthogonal invariant, is known as the Gaus- 
sian orthogonal ensemble (GOE). Likewise, the case, where the entries of 
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X are i.i.d. Gaussian, is known as the (null or white) Wishart ensemble. In 
particular, the Wishart ensemble has been used in statistics since the 30s as 
an important element of the sample covariance analysis, the principal com- 
ponent analysis first of all, in the asymptotic regime n oo,m < oo (see, 
e.g., [21] and the references therein). The eigenvalue distribution of these 
matrices for n — > oo, m — > oo, m/n c£ [0, oo), that is, an analog of the law 
of large numbers for n~^J\fn[^], was found in [20]. 

The CLT for certain linear eigenvalue statistics of the Wigner and the sam- 
ple covariance matrices were also considered in recent papers [2, 5] and [8]. 
In [2] the Wigner and the sample covariance matrices (in fact, more general 
matrices) and linear eigenvalues statistics for polynomial test functions were 
studied by using a considerable amount of nontrivial combinatorics, that is, 
in fact, a version of the moment method of proof of the CLT. This requires 
the existence of all moments of entries and certain conditions on their growth 
as their order tends to infinity. The conditions were then relaxed for differ- 
entiable test functions under the additional assumption that the probability 
law of entries satisfies a concentration inequality of the Poincare type. 

Related results are obtained in [8] for a special class of Wigner matrices, 
whose entries have the form 

Wjk = u{Wjk), ueC^R), 

(1-3) 

sup |n (x)! < oo, sup|ii (x)|<oo, 

where {Wjk}i<j<k<n are the independent standard (ElVTjfc} = 0, E{W"J^} = 
1) Gaussian random variables. For these, rather "close" to the Gaussian, 
random matrices a nice bound for the total variation distance between the 
laws of their linear eigenvalue statistics and the corresponding Gaussian 
random variable was given. The bound is then used to prove the CLT for 
linear eigenvalue statistics with entire test functions without explicit formula 
for the variance with a possibility of extension to functions and also for 
certain polynomials of growing with n degree, as in [29] . 

In [5] the real symmetric and Hermitian sample covariance matrices (in 
fact, more general matrices) were studied, assuming that the entries Xaj, 
a = 1, . . . ,m, j = 1, . . . ,n, of X are such that 

(1.4) E{Xn} = 0, B{Xf,} = l, E{Xfj = 3 
in the real symmetric case and 

(1.5) E{Xn} = E{x2j = 0, E{|Xn|2} = l, E{|Xn|^} = 2 

in the Hermitian case. Under these conditions the CLT for linear eigenvalue 
statistics with real analytic test functions was proved. 
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Conditions (1.4) and (1.5) mean that the fourth cumulant (known in 
statistics as the excess) of entries is zero. On the other hand, it was shown 
in [18] that for ^p{\) = (A — z)~^,'isz ^ the variance of the corresponding 
hnear statistic (the trace of resolvent) of Wigner matrices contains the fourth 
cumulant of entries. Thus, even in the class of real analytic test functions 
one can expect more in the case of nonzero fourth cumulant of entries. 

The requirement of the real analyticity of test functions results from the 
use of the Stieltjes transform of the eigenvalue counting measure as a ba- 
sic characteristic (moment generating) function. The Stieltjes transform was 
introduced in the random matrix studies in [20] and since then proved to be 
useful in a number of problems of the field (see, e.g., [4, 12, 18] and [23] and 
the references therein). We found, however, that while studying the CLT of 
the above ensembles it is more convenient to use as a basic characteristic 
function not the collection of moments or the Stieltjes transform but the 
Fourier transform of the eigenvalue counting measure, that is, the standard 
characteristic function of probability theory. This allows us to prove the CLT 
for linear eigenvalue statistics with sufficiently regular (essentially C^) test 
functions (but not real analytic as in [5]) and assuming the existence of the 
fourth moments of entries satisfying a Lindeberg type condition (but not all 
the moments or the Poincare type inequality as in [2], conditions (1.3) as 
in [8], or conditions (1.4) and (1.5) as in [5]). Besides, all proofs follow the 
same scheme based on the systematic use of rather simple means: the Fourier 
transform, certain differential formulas, that is, a version of integration by 
parts [see (2.20) and (3.6)], and an "interpolation trick" [25], which allows 
us to relate the asymptotic properties of a number of important quantities 
for general entries and those for the Gaussian entries. For both classes of 
random matrices we prove first the CLT for matrices with Gaussian entries 
(the GOE and the Wishart ensemble, see Theorems 2.2 and 4.2), then con- 
sider matrices with zero excess of entries, where the CLT can be obtained 
practically directly from that for the GOE and the Wishart ensemble by 
using an "interpolation" trick (see the proof of Theorems 3.3, 3.4 and 4.3). 
Finally, the proofs in the general case of nonzero excess of entries essentially 
follow those of the GOE and the Wishart cases and use again the interpola- 
tion trick that makes the proofs shorter and simpler (see Theorems 3.6 and 
4.5). 

The paper is organized as follows. In Section 2 we present the basics of our 
approach by proving the CLT for linear eigenvalue statistics in a technically 
simple case of the GOE (for other proofs see, e.g., [7, 12, 13] and [15] and 
the references therein). In Section 3 we consider the Wigner matrices and 
in Section 4 the sample covariance matrices. We confine ourselves to real 
symmetric matrices, although our results as well as the main ingredients of 
the proofs remain valid in the Hermitian case with natural modifications. 

Throughout the paper we write the integrals without limits for the inte- 
grals over the whole real axis. 
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2. Gaussian orthogonal ensemble. 

2.1. Generalities. We recall first several technical facts that will be often 
used below. We start from the generalized Fourier transform, in fact, the 7r/2 
rotated Laplace transform (see, e.g., [32], Sections 1.8-9). 

Proposition 2.1. Let f-.M.-^ be locally Lipshitzian and such that 
for some 6 > 

(2.1) supe-^^\f{t)\<oo 

t>o 

and let f :{z gC: '^z < —6} ^ C be its generalized Fourier transform 

(2.2) f{z) = i^' e-'^'f{t)dt. 



The inversion formula is given by 







(2.3) f^t) = —J^e^^'f{z)dz, t>0, 

where L = {—oo — i£,oo — ie), e > 5, and the principal value of the integral 
at infinity is used. 

Denote for the moment the correspondence between functions and their 
generalized Fourier transforms as f ^ f . Then we have the following: 

(i) /'(t)^i(/(+0) + z/(z)); 

(ii) f',f{T)dT^{iz)-^f{z)- 

(iii) /o* /i(t - r)/2(r) dr := (/i * f^m ^ ifi{z)f2{z); 

(iv) if P , Q and R are locally Lipshitzian, satisfy (2.1), and 

(2.4) l + iQ(z)/0, 9z<0, 
then the equation 

(2.5) P{t)+ I Q{t-ti)P{ti)dti = R{t), t>0, 

Jo 

has a unique locally Lipshitzian solution 

(2.6) P^R{l + iQ)-^ 
or 

(2.7) P{t) = -i f T{t-ti)R(ti)dti, 

Jo 

where 

(2.8) T^{l + iQ)-\ 
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In particular, if R{t) is differentiable, R{0) = 0, and 

(2.9) Q{t)= f Qi{s)ds, 

Jo 

then the equation 

(2.10) P{t)+ f dti r Qi(ti-t2)P{t2)dt2 = R{t), t>0, 

Jo Jo 

has a unique locally Lipshitzian solution 

(2.11) P{t) = - f Ti{t-ti)B!{ti)dtu 

Jo 

where 

(2.12) Ti^{z + Qi)-^ 
provided by 

(2.13) z + Qi{z)^0, 9z<0. 

The next proposition presents a simple fact of linear algebra: 

Proposition 2.2 (Duhamel formula). Let Mi, M2 be n x n matrices 
and t£M. Then we have 

(2.14) e(M^+M2)t ^ ^M^t ^ /•*^Mi(f-s)^^g(Mi+M2)s^^_ 

Jo 

Consider now a real symmetric matrix M = {Afjfc}^^^]^ and set 

(2.15) C/(^) = e^*^^ t£R. 

Then U{t) is a symmetric unitary matrix, possessing the properties 
Uit)U{s) = U{t + s), 

(2.16) 

\\U{t)\\ = l, \U,k{t)\<l, 5]|C/,fc(t)|2 = l. 

The Duhamel formula allows us to obtain the derivatives of U{t) with respect 
to the entries Mjk, j,k = 1, . . . ,n, of M: 

(2.17) D.kUabit) = if3jk[{Uaj * U,k){t) + {Ubj * f/afc)(t)], D^k = d/dMjk, 
where 

(2-18) /?,fc = (i + 5,.)-^ = {5/^' 
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and the symbol "*" is defined in Proposition 2.1(iii). Iterating (2.17) and 
using (2.16), we obtain the bound 

(2.19) \D]^Uab{t)\<ci\t\\ 

where ci is an absolute constant for every /. 

The next proposition presents certain facts on Gaussian random variables. 



Proposition 2.3. Let ^ = he independent Gaussian random 

variables of zero mean, and ^-.W he a differentiate function with poly- 
nomially hounded partial derivatives / = 1, . . . ,p. Then we have 

(2.20) E{6$(0} = mtmm)}^ i=i,...,p, 

and 

(2.21) Var{a>(e)}<EE{e?}E{|$KOl'}- 

1=1 

The first formula is a version of the integration by parts. The second is a 
version of the Poincare inequality (see, e.g., [6]). 

Next is the definition of the Gaussian orthogonal ensemble. This is a real 
symmetric n x n random matrix 

(2.22) M = n-^/^W, W = {W,k eR,Wjk = Wk,}lk=i, 
defined by the probability law 

(2.23) Z-/e-T^^V4«.^ Yl dW,,, 
where Z„i is the normalization constant. Since 

l<i<" i<j<k<n 

the above implies that {VFjfc}i<j<fc<„ are independent Gaussian random 
variables such that 

(2.24) B{Wjk} = 0, B{W^f,} = w\l + 6jk). 

Here is a useful bound for linear eigenvalue statistics of the GOE [9, 23]: 

Proposition 2.4. Let M he the GOE matrix (2.22)-(2.24) and Mn['^] 
he its linear eigenvalue statistic (1-1), corresponding to a differentiahle test 
function. Then 

(2.25) Var {A/; [(/?]} < 2w^'¥.{n^^ Ti {M){^' {M)*)] 

(2.26) <2w^Un^\ip" 
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Proof. The spectral theorem for real symmetric matrices and (1.1) 
imply that 

(2.27) MnW]=TT^{M). 

Thus, we can apply (2.21) to <1>(M) = TTip{M), viewing it as a differentiable 
function of the independent Gaussian random variables Mjk = n~^/'^Wjk, 
^ < 3 <k <n, satisfying (2.24). By using (2.21), (2.24) and the formula [see 
(2.17)] 

(2.28) TV v?(M) = 2/J,-fc(/p;.fe(M), 
we obtain 

Var{A4M} <u;2n-i E 4(1 + <5,-fc)/?|,E{|(^;.,(M)|2} 

i<i<fc<" 

n 

= 2w^n-^ E E{|v.;-,(M)|2} = 2u;2E{n-iTr^'(Af)(^'(M))*}. 

j,k=i 

This yields (2.25). Using it and the inequalities 

(2.29) |TrA| <np||, ||V'(5)|| < sup |^/'(A)|, 

AeM 

valid for any normal matrix A, Hermitian matrix B, and -i/^rR— >C, we 
obtain (2.26). □ 



We recall now an analog of the law of large numbers for linear eigenvalue 
statistics of the GOE (see, e.g., [4, 12] and [23] and the references therein). 

Theorem 2.1. Let M be the GOE matrix (2.22)-(2.24), and Mn[(p] be 
a linear statistic of its eigenvalues (1.1). Then we have for any bounded and 
continuous : M — > C with probability 1 

(2.30) \mi^n-^MnM= J v{X)N,^i{d\), 
where the measure 

(2.31) N,,i{dX) = psci(A) dX, psci(A) = (2W)-i(4^2 _ ^2)i/2 
is known as the Wigner or the semicircle law, and = max{0,a;}. 

We need below a weak version of the theorem in which the convergence 
with probability 1 is replaced by the convergence in mean, that is, 

(2.32) hm E{n-^MrM}= I v{X)N,,,{dX) 
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for any continuous and bounded (p. We outline the proof of this relation 
to introduce several elements of techniques used in the paper (see [23] for 
details). 

Introduce the normalized counting (empirical) measure of eigenvalues as 

(2.33) iV„(A) = KaS"^ G a : / = 1, . . . , n]/n. 

Then we have 

E{n-iA/:„M}= / (^(A)E{iV„(dA)}, 



hence, (2.32) is equivalent to the weak convergence of ^{N^} to N^^i. More- 
over, since by (2.24) 

(2.34) j \^^{Nn{dX)] = ¥.{n-^TTW'^] = {l + n^^)w'^, 

the sequence {EjA^^^}} is tight, and it suffices to prove the vague convergence 
of the sequence, for instance, the convergence of the Stieltjes transforms 

p.35) = /5M 

of EjA'n} for any SJz / (see, e.g., [1], Section 59) to the Stieltjes transform 
(2.36) f{z) = (V^^ _ 4y;2 _ ^) /2^2 



of A'sci, where \J — Aufl is defined by the asymptotic \J — 4w'^ = z + 
0{z^^), z ^ oo. 

It follows from the definition of Nn that 

(2.37) /„(z)=E{n-iTrG(z)}, 
where 

G{z) = {M-z)-\ Qz^O, 
is the resolvent of M. We will need the resolvent identity 

(2.38) {A + B- z)-^ -{A- z)-^ = -{A + B- zy^B{A - zy\ 



its implication 

d_ 

de 

and the bounds 



(2.39) -j-^{A + eB- z)^^ = -{A - z)-^B{A - z)-^ 



(2.40) ||(^ - zY'W < \^z\-\ \{{A - z)-')^,\ < \-sz\-\ 
valid for real symmetric matrices A and B. 
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It follows from (2.38) for A = 0,B = M that 



^<j<k<n 



where Pjk is defined in (2.18). Since Wjk,l < j < k < n, are independent 
Gaussian variables, we can write, in view of (2.20) and (2.24), 



1 2?/;^ r 1 
(2.41) ^"^^) = -i + ^^| E D,,G,,iz)j, 



1 2u;2. 

where Djk is defined in (2.17). It follows from (2.39) that [of. (2.17)] 

(2.42) DjkGabiz) = -Pjk{Gaj(z)Gkb{z) + Gak{z)Gjb{z)). 
This allows us to write (2.41) as 

(2.43) Uz) = -Z-' - w'z-'B{gl{z)} - w^z~'B{7^-^ Tr G\z)}, 
where 

(2.44) gn{z)=n-^TrGiz). 

By using (2.26) with 93(A) = (A — z)~^, we find that 

2w^ 



(2-45) Var{g4z)}< 



hence, 



|E{g^W}-/n(^)l<Var{5„(z)}< 



Besides, (2.29) and (2.40) imply that \'E{n~^TrG^{z)}\ < l/n\Qz\^, and 
(2.35) implies that |/„| < l/\Qz\. 

In view of the above bounds, the sequence {/«} is compact with respect 
to the uniform convergence on any compact set K cC\M. and the uniform 
on K limit / of any convergent subsequence of {/„} satisfies the quadratic 
equation 

(2.46) fiz) = -z-'-w^z-'fiz), 

following from (2.43). In addition, we have, by (2.35), Qfn{z)Qz > 0, thus, 
Qf{z)Qz > 0. Now it is elementary to check that the unique solution of the 
above quadratic equation that satisfies this condition is (2.36). 
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2.2. Central limit theorem for linear eigenvalue statistics of differentiable 
test functions. The CLT for the GOE was proved in a number of works 
(see [7, 13] and [23] and the references therein). We present below a proof, 
whose strategy dates back to [18] and is used in the proofs of other CLT of 
the paper. 



Theorem 2.2. Let (p:M.^M. be a bounded function with bounded deriva- 
tive, and Mni^p] be the corresponding linear eigenvalue statistic (1.1) of the 
GOE (2.22)~(2.24). Then the random variable 

(2.47) M°[^]=Mn[^]-B{Mn[ip]} 

converges in distribution to the Gaussian random variable with zero mean 
and variance 

where 

(2.49) Aip = ip{Xi)-^{X2), AA = Ai-A2. 



Proof. By the continuity theorem for characteristic functions, it suf- 
fices to show that if 

(2.50) Z„(2;) = E{e^^^"[^l}, 
then for any x G M 

(2.51) Um ZJx) = Z(x), 

n— ♦CO 

where 

(2.52) Z{x) = exp{-xVGOE[^]/2}. 

We obtain first (2.52), hence the theorem, for a certain class of test functions 
and then we extend the theorem to the bounded functions with bounded 
derivative, by using a standard approximation procedure [see also Remark 
2.1(2) for a wider class]. 

Assume then that ip admits the Fourier transform 

(2.53) ^^^^ = hj ^"'^"PWdX, 
satisfying the condition 

(2.54) /(l + |t|2)|<^(t)|dt<cx). 
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Following the idea of [18], we obtain (2.52) by deriving the equation 

(2.55) Z{x) = l- Fgoe M r yZ{y) dy. 

Jo 

The equation is uniquely soluble in the class of bounded continuous functions 
and its solution is evidently (2.52). 
It follows from (2.50) that 

(2.56) Z;(rr) = iE{AA°[(^]e"^"[^]}. 
This, the Schwarz inequality and (2.26) yield 

\Z'^{x)\<V2w(snp\^'{X)\). 
Since Zn{0) = 1, we have the equality 

(2.57) Zr,{x) = l+ r Z'^{y)dy, 

Jo 

showing that it suffices to prove that any converging subsequences {Zn-} 
and {Z'n^} satisfy 

(2.58) lim ZnAx) = Z{x), lim Z' {x) = -xVgof.Z{x). 

Indeed, if yes, then (2.58), (2.57) and the dominated convergence theorem 
imply (2.55), hence (2.52). 

The Fourier inversion formula 

(2.59) ^(A) = j e'*^(p{t) dt 
and (2.27) yield for (2.47) 

(2.60) M°M= j m<{t)dt, 
where 

(2.61) Un{t) = TiU{t), <(t)=n„(t)-EK(t)}, 

and U{t) defined by (2.15) with the GOE matrix M. Thus, we have by (2.56) 

(2.62) Z'^{x)=i [ 0{t)Yn{x,t)dt, 



where 

(2.63) Yn{x,t) = EK(t)e„(x)}, e„(x) = e^^-^-M. 

Since 



(2.64) Yn{x,t)=Yn{-x,-t) 
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we can confine ourselves to the half-plane {t>0,x£ M}. 

We will prove that the sequence {Yn} is bounded and equicontinuous on 
any finite set of {t >0,x € ffi}, and that its every uniformly converging on 
the set subsequence has the same limit Y, leading to (2.58). This proves 
the assertion of the theorem under condition (2.54). Indeed, let {Zni}i>i be 
subsequence converging to Z Z. Consider the corresponding subsequence 
{Yni}i>i- It contains a uniformly converging sub-subsequence, whose limit 
is again Y, and this forces the corresponding subsequence of {Zni}i>i to 
converge to Z, a contradiction. 

It follows from (2.25) and (2.26) with ip{X) = e**^ and ip{X) = iXe'*^ and 
from (2.24) that 

(2.65) Var{unit)} <2w'^t'^ 
and 

(2.66) Var{<(t)} < 2w'^n^^E{Tiil +t^M^)} < 2w'^{l + 2w'^t'^). 
This, (2.63), the Schwarz inequality, and |en(x)| < 1 yield 

(2.67) \Yn{x,t)\ < E{|n°(t)|} < V2w\t\ 
and 

d 



(2.68) 



dt 



Yn{x,t) 



< Vari/2{n;(t)} < V2w{l + 2^1-2*2)1/2 
and according to the Schwarz inequality, (2.26) and (2.65), 

= |EK(t)AA°[(^]e„,(x)}| 

< Var^/2|^^(^)| Yari/2{_A/;[^]} <2w'^\t\ sup|(^'(A)|. 



(2.69) 



We conclude that the sequence {Yn} is bounded and equicontinuous on any 
finite set of M^. We wih prove now that any uniformly converging subse- 
quence of {Yn} has the same limit Y, leading to (2.55), hence to (2.51) and 
(2.52). 

It follows from the Duhamel formula (2.14) with Mi = and M2 = iM 
and (2.61) that 



Un{t)=n + i / V MjkUjkih)dti, 



hence, 
(2.70) 



Yn{x,t) 



/ ^{WjkU,k{ti)e°n{x)}dh, 
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where e° = e„ — E{e„} or, after applying (2.20) and (2.24), 

2 /-t " 

Now, by using (2.17) and (2.28), we obtain that [cf. (2.42)] 



(2.71) y„,(a;,t) = / ^ (1 + <^,fc)E{I),fe([/,fc(ti)e° (x))} dii 

■^o ,-,fc=l 



(2.72) L>jA..e„(x) = 2i/3jfexe„(x)v?jfc(M) = -2/?jfcxen(x) / tUjk{t)(p{t) dt, 
where the last equality follows from [see (2.59)] 

(2.73) if'iM) = i f (f{t)tU{t) dt. 



This and (2.71) yield [cf. (2.43)] 

Yn{x,t) = -w'^n-^ f tiE{n„(ti)e°(rE)}dti 
Jo 



uP'n ^ 



dti [ E{Un{t2 - ti)Un{t2)e°Xx)} dt2 

Jo 

- 1v?x r E{e„(x)n^^ Tr U{txy(M)\ dti. 
Jo 

Writing 

(2.74) Vn{t)=n~'B{unit)} 
and 

(2.75) Un{t) =u°^{t) +nvn{t), e„(x) = e°(x) + Z„(x), 
we present the above relation for Yn. as 

Ynix,t) + 2w'^ f dti r Vn{tl-t2)Yn{x,t2)dt2 

(2.76) ^° ^° 

= xZn{x)An{t) +rn{x,t) 

with 

(2.77) An{t) = -2w;2 /* E{n~^ Tr U{ti)ip'{M)] dti 

Jo 

and 

rn{x,t) = -v?n^^ I tiYn{x,ti)dti 
Jo 

(2.78) - w'^n-^ r dti r E{<(ii - t2)<(i2)e° (x)} dt2 

Jo Jo 

-liw'^xn-^ dti J t2'f{t2)B{un{ti + t2)el{x)} dt2 
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where we used (2.16) and (2.73). 

It follows from the inequahty |e°(x)| < 2, the Schwarz inequahty, (2.54) 
and (2.65) that the hmit 

(2.79) lim rJx,t) = 

n — ^oo 

holds uniformly on any compact of {t>0,x£ M}. Besides, by Theorem 2.1 
the sequences {vn} of (2.74) and {An} of (2.77) converge uniformly on any 
finite interval of M to 

1 i'2w 

(2.80) vit) = ^ / e'*V4u;2_ A2(iA 

znw^ J-2w 

and 

1 rt p2w 

(2.81) A{t) = -- / dh / e**iV(A)V4Ti;2-A2dA. 

vr Jo J-2w 

The above allows us to pass to the limit n; ^ oo in (2.76) and obtain that 
the limit Y of every uniformly converging subsequence {1^, }/>i satisfies the 
equation 

(2.82) Y{x,t) + 2w^ f dti r v{ti-t2)Y{x,t2)dt2 = xZ{x)A{t). 

Jo Jo 

The equation is of the form (2.10), corresponding to 5 = in (2.1), thus, we 
can use formulas (2.12) and (2.13) to write its solution. 

It follows from (2.2) and (2.80) (or from the spectral theorem and Theo- 
rem 2.1) that 

(2.83) v = f, 

where / is the Stieltjes transform (2.36) of the semicircle law. Thus, in our 
case (2.4) takes form 



(2.84) z + 2w^f{z) = Vz^ - Aw^ / 0, '^z^O, 

and 



^'^'^ = ^JLZ + 2w^fiz)''- 

Replacing here the integral over L by the integral over the cut [—2w, 2w] 
and taking into account that \/ z^ — A.nP' is ibi\/4t(;2 — A2 on the upper and 
lower edges of the cut, we find that 
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Besides, in our case the r.h.s. of (2.10) is xZ{z)A{t), hence, we have by (2.11) 
ixZ(x) f^^ da 



vr^ j-2w JAw"^ - u? 
(2.86) 

^- — ^ip'{X)VAw^ - X'^dX, t>0. 

-2w (A — /i) 

According to (2.64), the same hmiting expression is valid for i < and x G '. 
thus, we have, in view of (2.54) and (2.62), 

xZ{x) f"^^ dfi 



2w \/ AvS^ — 



2w 



Writing 



^^^M^^ ^^(A)V4^^^dA. 

2w (A-/U) 



15 2 



y,'{X){^{X) - (^(^)) = 2^(V^(A) - 



and integrating by parts with respect to A, we obtain 

f2 87) lim Z' (x)- r n^^y" iiw'-Xf,)dXdt, 

hence (2.58), and then (2.55), thus the assertion of the theorem under the 
condition (2.54). 

The case of bounded test functions with bounded derivative can be ob- 
tained via a standard approximation procedure. Indeed, for any bounded <p 
with bounded derivative there exists a sequence {(fk} of functions, satisfying 
(2.54) and such that 

sup|v3fe(A)| <sup|(/?'(A)|, 

AgK AeK 

(2.88) 

lim sup |v?'(A) -(/9fc(A)| =0 > 0. 

By the above proof we have the central limit theorem for every ip^- Denote 
for the moment the characteristic functions of (2.50) and (2.52) by Zn[f] 
and Z[cp], to make explicit their dependence on the test function. We have 
then for any bounded test function with bounded derivative 

\Zn[ip] - Z[ip]\< \Zn[(p] - Zn[fk] \ + \Zn[Pk] " Z[ipk]\ 

(2.89) 
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where the second term of the r.h.s. vanishes after the hmit n — > oo in view 
of the above proof. It follows from (2.25), (2.33) and (2.50) that 

< V2w\x\ l^j - ^UA)pE{iV„(dA)} 

Now (2.88) implies that the integral on the r.h.s. is bounded by 

2sup|(^'(A)|2E{iV„(M\[-AA])}+ sup y{\)-^'k{X)\\ A>2w, 
Asm |A|<yl 

where the first term vanishes as n — > oo by (2.32), and the second term 
vanishes as A: — > oo by (2.88). Besides, according to (2.52), 

|ZM - Z[^k]\ < x^\VgoeM - VGOE[^k]\/2 

and taking into account the continuity of Vgoe of (2.48) with respect to the 
convergence on any interval |A| < A, A> 2w, we find that the third term 
of (2.89) vanishes after the limit k oo. Thus, we have proved the central 
limit theorem for bounded test functions with bounded derivative. For wider 
classes of test functions see [7] and [15] and Remark 2.1. □ 

Remark 2.1. (1) We note that the proof of Theorem 2.2 can be easily 
modified to prove an analogous assertion for the Gaussian unitary ensemble 
of Hermitian matrices, defined by (2.22) with Wjk = W^j and the probability 
law 

n n 
j=l i<j<k<n 

The result is given by Theorem 2.2 in which Vqoe is replaced by Vgue = 
Vgoe/2. 

(2) It follows from the representation of the density p„ of ^{Nn} via the 
Hermite polynomials (see [19], Chapters 6 and 7) or from Theorem 2.3 of 
[26] that 

for finite c > 0, C < oo, and a sufficiently big |A|. This bound and the ap- 
proximation procedure of the end of proof of Theorem 2.2 allows one to 
extend the theorem to test functions whose derivative grows as Cie'^^^ 
for any ci > and Ci < oo. 
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3. Wigner ensembles. In this section we prove the central limit theorem 
for linear eigenvalue statistics of the Wigner random matrices. We start from 
the analog of the law of large numbers, by proving that the normalized count- 
ing measure of eigenvalues converges in mean to the semicircle law. The fact 
has been well known since the early seventies (see [4, 12] and [22] and the 
references therein). We give a new proof, valid under rather general condi- 
tions of the Lindeberg type and based on a certain "interpolation" trick that 
is systematically used in what follows. We then pass to the CLT, proving it 
first for the Wigner ensembles, assuming the existence of the fourth moment 
of entries, their zero excess and the integr ability of (1 -|- |i|^)<^ (Theorems 
3.3 and 3.4), and then proving it in the general case of an arbitrary excess 
(Theorem 3.6), assuming the existence of the fourth moments satisfying a 
Lindeberg type condition (3.57) and again the integrability of (1 -|- |i|^)(^. 

3.1. Generalities. We present here the definition of the Wigner ensem- 
bles and technical means that we are going to use in addition to those given 
in the previous section. 

Wigner ensembles for real symmetric matrices can be defined as follows: 

(3.1) M = n-^/^W, W = {H^j,") G M, W^^^ = t^if 

[cf. (2.22)], where the random variables Wj'^\ I < j < k < n, are indepen- 
dent, and 

(3.2) B{W^^^} = 0, E{(W,5))2} = (1 + 5,,.)^^ 

that is, the two first moments of the entries coincide with those of the GOE 
[see (2.22)-(2.24)]. In other words, the probability law of the matrix W is 

(3.3) P{dW)= n Fjk\dW,,), 

^<3<k<n 
In) 

where for any 1 < j < k < n, F\ is a probability measure on the real line, 
satisfying conditions (3.2). We do not assume in general that is n- 

independent, and that F^-^ are the same for 1 < j < k <n and for j = k = 
1, . . . ,n, that is, for off-diagonal and diagonal entries as in the GOE case. 

Since we are going to use the scheme of the proof of Theorem 2.2 (the 
CLT for the GOE), we need an analog of the differentiation formula (2.20) 
for non-Gaussian random variables. To this end, we recall first that if a 
random variable ^ has a finite pth. absolute moment, p >1, then we have 
the expansions 

/(t):=E{e^*n = E^(^tF+0(O 
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and 

l{t) := logE{e**«} = f: ^iity + oin, t ^ 0. 

Here "log" denotes the principal branch of logarithm. The coefficients in 
the expansion of / are the moments {/ij} of ^, and the coefficients in the 
expansion of / are the cumulants {kj } of ^. For small j one easily expresses 
Hj via //I, //2, • • • , Aij- In particular, 

^^^^ Kl=/ii, K2 = ^2 -/^? = Var{^}, K3 =/i3 - 3/X2/il + 2;uf, 

K4 = /i4 - - Afl-sHi + l2fJ,2^4 - 6/ii, .... 

In general, 

(3.5) Kj = ^cxn\, 

A 

where the sum is over all additive partitions A of the set {1, . . . , j}, cx are 
known coefficients and fix = HieA^i; see, for example, [27]. 
We have [18]: 

Proposition 3.1. Let ^ he a random variable such that E{|^|*'+^} < oo 
for a certain nonnegative integer p. Then for any function $ : M — > C of the 
class C^^^ with bounded derivatives Z = 1, . . . ,p + 1, we have 

(3.6) nmo} = E ^E{$W(e)} + e„ 

1=0 

where the remainder term £p admits the bound 

(3.7) \s,\<C,E{\Cr'}snp\^^^+'\t)l < ' ^ f ! ^ • 

Proof. Expanding the left- and the right-hand side of the identity 

B{^e^'i} = f{t)nt) 

in powers of it, we obtain 

(3.8) iJ.r+i = ^r.]nj+ifJ.r-j, r = 0,l,...,p. 

j=o ^-^ ^ 

Let vr be a polynomial of degree less or equal p. Then (3.8) implies that (3.6) 
is exact for $ = vr, that is, is valid with ep = 0: 



E{MO}=E^E{vr(^'^(0}- 



j=o -J 
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In the general case we write by Taylor's theorem $ = vr^ + r^, where vTp is 
a polynomial of degree p and 

rp(t) = ^ [\^+^\tv){l-vYdv. 



p! Jo 



Thus, 

Kq> 
(p+1)! 

where 



(3.9) - m^m < mrpim < jz-^Bu^r'}, 



i^$=sup|$(P+^)(t)| <CX). 



Besides, 



(3.10) 



^(^\t) - Tr'ij^t) = f^^,[, f\'^P+^\tv){l-vy-^dv, l = 0,...,p, 



and, therefore 



E{evrp(0}-E^E{$0)(e)} 



?! 



The sum on the r.h.s. can be estimated with the help of the bound [27]: 

(3.11) \K,\<fE{\C-E{CW}. 

Since (a + by < 2^~^{a^ + b^) for a positive integer j and nonnegative a and 
6, we have 

(3.12) < i^E{(iei + imm < m^nm- 

This bound and the Holder inequality E{|^p} < E{|^|p+2}jV(p+2) yjeld 
.t^o i!(p-J + l)! - ^' ' ^,eoj!(p-J + l)! 

(3.13) 

- ^ (p+l)! 

The proposition now follows from (3.9)-(3.13). □ 



Here is a simple "interpolation" corollary showing the mechanism of prox- 
imity of expectations with respect to the probability law of an arbitrary 
random variable and the Gaussian random variable with the same first and 
second moments. Its multivariate version will be often used below. 
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Corollary 3.1. Let ^ be a random variable such that E^{|^|p+^} < co 
for a certain integer p> 1, E^j,^} = 0, and let be the Gaussian random 
variable, whose first and second moments coincide with those of ^. Then for 
any function ^> : R ^ C of the class C'P'^'^ with bounded derivatives, we have 

^dm} - ^^m} = E ITT / ^{$('+1) m)}s(^-'y' ds + 4, 

(3.14) 

where the symbols E^{...} and E^...} denote the expectation with respect 

to the probability law of ^ and^, {nj} are the cumulants of ^, E{. . .} denotes 
the expectation with respect to the product of probability laws of ^ and ^, 

(3.15) ^(s) = sV2^ + (l-s)V2^, 0<s<l, 

(3.16) |e;|<CpE{|^r2}sup|cI,(P+2)(t)| 

teK 

and Cp satisfies (3.7). 

Proof. It suffices to write 

(3.17) = f'^Bmmds 

Jo ds 

= \[' E{S-V2^CI>'(^(,)) - (1 - s)-l/2^$'(^(s))} ds 

^ Jo 

and use (3.6) for tlie first term in tlie parentlieses and (2.20) for the second 
term. □ 

3.2. Limiting normalized counting measure of eigenvalues. We will also 
need an analog of Theorem 2.1 on the limiting expectation of linear eigen- 
value statistics of Wigner matrices. It has been known since the late 50s that 
the measure is again the semicircle law (2.31) (see [4, 12, 22] and [23] for 
results and references). We give below a new proof of this fact [25] that is 
based on the matrix analog of the "interpolation trick" (3.14) and illustrates 
the mechanism of coincidence of certain asymptotic results for Gaussian and 
non-Gaussian random matrices. The trick will be systematically used in what 
follows. 

Theorem 3.1. Let M = n~'^/'^W be the Wigner matrix (3.1)-(3.3), sat- 
isfying the condition 

(3.18) u;3:=sup max E{|W^"^|^} < oo 

n i<j<k<n •' 
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and Nfi be the normalized counting measure of its eigenvalues (2.33). Then 

lim E{iV„} = iVsei, 

n— too 

where Nsd is the semicircle law (2.31) and the convergence is understood as 
the weak convergence of measures. 

Proof. It follows from (3.2) that we have (2.34) for the Wigner ma- 
trices. Thus, the sequence {E{A'^„}}„>o is tight, and it suffices to prove its 
vague convergence, or, in view of the one-to-one correspondence between 
the nonnegative measures and their Stieltjes transforms (see, e.g., [1]), it 
suffices to prove the convergence of the Stieltjes transform of expectation of 
the normalized counting measure [see (2.35) and (2.37)] on a compact set of 
C\M. Let M = n-^/^W be the GOE matrix (2.22)-(2.24), and G{z) be its 
resolvent. Then, by Theorem 2.1, it suffices to prove that the limit 

(3.19) lim |E{n-iTrG(z)} -E{n~iTrG(z)}| =0 

n — too 

holds uniformly on a compact set of C \ R. 

Following the idea of Corollary 3.1, consider the "interpolating" random 
matrix [cf. (3.15)] 

(3.20) M(s) = s^/2^/ + (1-s)^/2m, 0<s<1, 

viewed as defined on the product of the probability spaces of matrices W and 
W. In other words, we assume that W and W in (3.20) are independent. We 
denote again by E{. . .} the corresponding expectation in the product space. 
It is evident that M(l) = M, M(0) = M. Hence, if G{s,z) is the resolvent 
of M(s), then we have 

n-^E{TVG(z) -TVG(z)} 

(3.21) = -^'E{n'^TrG{s,z)}ds 

Jo OS 

where we used (2.39) and (3.20). 

Now we apply the differentiation formula (3.6) to transform the contri- 
bution of the first term in the parentheses of the r.h.s. of (3.21). To this 
end, we use the symmetry of the matrix {Gjk} to write the corresponding 
expression as 

(3.22) (n3.)-V2 ^ (3^,B{W^^\g%,}, 
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where f3jk are defined in (2.18) and we denote here and below 

G' = |g(„). 

(n) 

Since the random variables Wjf, , 1 < j < k < n, are independent, we can 
apply (3.6) with p = 1 and <J> = G^-^. to every term of the sum of (3.22). We 
obtain, in view of (3.2), (3.18) and (3.20), 

2 

(3.23) ^ J2 E{I),fc(s)(G')fcj}+ei, Dj^is) = d/dUj^is), 

i<i<fc<" 

where [cf. (3.7)] 

(3.24) kil<^ E sup \D%{s){G%,\, 

Sn is the set of n x n real symmetric matrices, and Ci is given by (3.7) for 
p = l. 

On the other hand, applying to the second term in the parentheses of 
(3.21) the Gaussian differential formula (2.20), we obtain again the first 
term of (3.23). Thus, the integrand of the r.h.s. of (3.21) is equal to ei. 

It follows from (2.42) and its iterations that 

(3.25) \D',kG,k\ < ci/\Qzf+'\ \D'^kiG'U\ < ci/\Qz\^'+^\ 

where q is an absolute constant for every /. The bounds and (3.24) imply 

CW3 



£1 < 



and C denotes here and below a quantity that does not depend on j, k and 
n, and may be distinct on different occasions. 

This and the interpolation property (3.21) yield the assertion of the the- 
orem. □ 

In fact, we have more (see [4, 12] and [22] for other proofs and references). 



Theorem 3.2. The assertion of Theorem 3.1 remains true under the 
condition 

n . 

(3.26) lim n-^ Y / F^':!\dW) = Vr > 0. 



j,k- 
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Proof. Given r > 0, introduce the truncated matrix 



(n)T 



(3.27) 

Let ii] il^ljk) and k^^^, {i^l.jk) be the ^th moment and cumulant of W^'^ 
(wj^^), respectively. Then 

W\'F'^l\dW). 



\f4,jk - ^J'l,jk\ < 2 / 

This and (3.5) yield 
(3.28) 



\W\>T^ 



ljk-'^iJk\<C f \W\^Fj'^\dW), 



l\W\>T^ 

where C depends only on Z. In particular, we have 



(3.29) 

and 



C 



T\ n 



\W\>T^ 



3 



(3.30) < 



C 



4~l 



W^F^'^\dW), l<4. 



As in the previous theorem, it suffices to prove the limiting relation (3.19). 
To this end, we show first that, for every r > 0, the limit 



(3.31) 



lim |E{ri-^TrG(z)} - E{n~^ Tr C {z)}\ =0 



with G'^{z) = {A'F^ — z) ^ uniform on any compact set of C \ M. Indeed, 
we have by the resolvent identity (2.38) and the bound \{G'^{z)G{z))ji:\ < 



\E{n-^Tr{G{z)-G^{z))}\ 



1 

-i^ y: b{{g^z)g{z)),,{w^:^-w^;:^^)} 



j,k=i 



< 



1 



77,3/2 1 3=2; p 



\W\Fj]:\dW)< 



where 
(3.32) 



n „ 

L„(T)=n-2 y: / 



j,k=l 



l\W\>T^ 



W'Fj-\dW). 



The last inequality and (3.26) imply (3.31). 
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Hence, it suffices to show that 

(3.33) Rnr ■■= E{?i~^ Tr - Ejn"^ Tr G(z)} 
vanishes after the subsequent hmits 

(3.34) n^oo, r^O, 

uniformly on any compact set of C \ M. 
Introduce the interpolation matrix 

(3.35) M''{s) = s^/'^M'' + sf/'^M, 0<s<l 

[cf. (3.20)], denote its resolvent by G(s, z) = {M'^{s) — z)~^ , and get an analog 
of (3.21): 

1 /•! " 

Rnr = -^, / E n{G'),u{s-'/'wlf^-{l-s)-'I^W,,)}ds. 

As in the previous theorem, we apply the differentiation formula (3.6) with 
p = 1 to every term containing the factors w'-'^'^ and Gaussian differentiation 
formula (2.20) to every term containing Wjk, and obtain 



Rr, 



1 



+ 4,A:E{D,fc(s)(G'),fc} 



j,k=l 



+ ^1 - — E ^'(1 + 5,k)^{Djk{s){G')jk} ds, 



where [cf. (3.24)] 



£1 < 



E n\w\tt} sup \D%{s){G'),,\< 



l<j<k<n M{s)eS„ 



and we took into account (3.25) and the bound 

B{\wj;:^^f}<Tn'/W{l + 6,,). 
Besides, we have by (3.25), (3.29) and (3.28) with / = 2, and the equalities 

Hl,jk = 0, K2Jk = W^(l + 6jk), 



n 



3/2 



E 4,,fcE{(G'),a 



j,k=i 



< 



C 



E (4jfc - ^'(1 + d,kmD,kis){G'),k} 



j,k=i 



< 



C 



Ln{T). 
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The last inequalities show that Rm- of (3.33) vanishes after the subsequent 
limits (3.34). This completes the proof of the theorem. □ 

Remark 3.1. Condition (3.26) is a matrix analog of the well-known 
Lindeberg condition 

in. 

(3.36) lim -V / x^Fl''\dx) = Vr > 

n^oonf-'J\x\>T,/E 



for a collection of independent random variables with probability 

laws Fj , j = l,...,n, to satisfy the central limit theorem. According to 
Theorem 3.2, the matrix analog (3.26) of the Lindeberg condition is sufficient 
for the validity of the semi-circle law for the Wigner ensembles. Thus, we can 
say that the semi-circle law is a universal limiting eigenvalues distribution 
of the Wigner ensembles in the same sense as the Gaussian distribution 
(the normal law) is universal for properly normalized sums of independent 
random variables. We mention two sufficient conditions for (3.26) to be valid, 
analogous to those of probability theory. The first is 

sup max E{|w|"^|2+^}<cx) 

for some 5 > 0. This is an analog of the Lyapunov condition of probabil- 
ity theory. The second sufficient condition requires that {l^jfc^}i<j<fc<n and 

^^^^^^^<j<'^ collections of independent identically distributed ran- 

dom variables, whose probability laws Fi and F2 do not depend on n and 
satisfy (3.2). This case generalizes the GOE, where Fi^2 are both Gaussian. 

3.3. Central limit theorem for linear eigenvalue statistics in the case of 
zero excess. We consider here a particular case of the Wigner ensembles, 
for which the fourth cumulant of entries, known in statistics as the excess, 
is zero. The case is of interest because here the limiting Gaussian law has 
the same variance as in the GOE case (Theorem 2.2), moreover, it can be 
obtained from that for the GOE by applying the "interpolation" trick that 
was used in the proof of Theorems 3.1 and 3.2 [see formulas (3.20) and 
(3.35)]. We start from an analog of Theorem 3.1. 

Theorem 3.3. Let M = n~^^'^W he the real symmetric Wigner matrix 
(3.1)-(3.3). Assume the following: 



(i) 



(3.37) w;5:=sup max ^{\w\l^\°} < oo- 
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(ii) the third and fourth moments do not depend on j , k and n: 

(3.38) ^3 = mwj'^^f}, /i4 = E{(W^5V}; 

(iii) the fourth cumulant of off-diagonal entries is zero: 

(3.39) K4 = /i4 — 3w^ = 0. 

Let (^iM^M be a test function whose Fourier transform (2.53) satisfies 
the condition 

(3.40) j {l + \t\^)\!p{t)\dt<oo. 

Then the corresponding centered linear eigenvalue statistic M^i^p] [see (2.47)] 
converges in distribution to the Gaussian random variable of zero mean and 
variance Vgoe[^] of (2.48). 

Remark 3.2. It may seem not too natural to have the {j,k) dependent 
second moments (3.2) of wj^'' and the (j, k) independent fourth moments 
(3.38). This is only for the sake of technical simplicity of the proof. In fact, 
it can be shown that the result does not depend on the diagonal entries, 
in particular, we can assume that the second moments will be the same for 
all 1 < j < k < n, or that only the fourth moments of off-diagonal entries 
are the same and the fourth moments of diagonal entries are just uniformly 
bounded. Likewise, we can replace (3.39) by 

lim sup max |«;4,fc|=0, 

n— >oo n l<j<k<n 

(n) 

where ^4^^ is the fourth cumulants of Wjf^ . 

Proof of Theorem 3.3. Let M = n'^l'^W be the GOE matrix (2.22)- 
(2.24) with the same variance of entries as the Wigner matrix, and MnYp\ 
be the centered linear eigenvalue statistic of the GOE. Then, in view of 
Theorem 2.2, it suffices to show that, for every x G M, 

(3.41) Rn{x) ■= E{e^^-^"[^]} - E{e^^-^"[^1} ^0, n ^ 00. 
Denoting 

(3.42) e„,(s,x) =exp{ixTrv3(M(s))°}, 
where M(s) is the interpolating matrix (3.20), we have 



Rn{x) = J^ — E{e„(s,x)}(is 



IX 



E{e° (s, rr) TV(^'(M(s))(5"i/2iy - (1 - .sy^/'^W)]ds 

Jo 
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[cf. (3.17)] or, after using (2.73), 

(3.43) Rn{x) = -^j\s J t^it)[An-Bn]dt, 



where 



1 " 



j,k=i 

1 " ^ 

(3.45) i?„ = E E{t^,.$4 

^) j,k=i 

with 

(3.46) ^'„ = e°(s,x)C/jfc(s,t), C/(s,t) = e^*^(^). 

Applying (3.6) with p = 3 to every term of (3.44) and (2.20) to every term 
of (3.45), we obtain (cf. Corollary 3.1) 

(3.47) An-Bn = T2 + T3 + es, 
where [cf. (3.23)] 

(3.48) Ti = ^fy^!y2 E ^i+i,jkE{D',kmn}, Djkis) = d/dM,kis), 

j,k=l 

and by (3.7) and (3.37) [cf. (3.24)], 

C' w " 

(3.49) 1^31-^72" XI \D%{s)^n\M{s)=M\ 



It follows then from (2.17), (2.19), (2.72) and (3.40) that 

(3.50) |L'jfc(s)$„| <Q(t,rE), 0</<4, 

and we denote here and below Ci{t, x) a polynomial in |t| and |x| of degree I, 
independent of j, k and n and not necessarily the same at each occurrence. 
This implies that 

(3.51) [esl <C4(t,x)/n^/2. 

Furthermore, it follows from (3.2) and (3.39) that n^jk = —96jkw'^. This, 
(3.48) for / = 3, and (3.50) yield 

(3.52) \T3\<C3it,x)/n. 
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To get a vanishing bound for T2 of (3.48), we use (2.17) and (2.72) to find 



the second derivative and we take into account (3.38) to have 



(3.53) 



■^T? E /3|feE|e°[(C/,fc*C/,fc*C/,fc)(t)+3(C/,fc*t/,,*t/fcfc)(t)] 
^ j,k=i ^ 

+ 2xen[{Ujk * Ujk){t) + {Ujj * Ukk){t)] 

X j e(p{e)Ujk{o)de 

2 



2x'enUjkit)i^J e(p{e)Ujkio)de^ 



where we write en for en(s,x), U{t) for U{s,t) and take into account that 
the convolution operation "*" of Proposition 2.1(iii) is commutative. 
Consider the two types of the sums above: 

n 

T2i=n-3/2 J2 Ujk{ti)Uyj{t2)Uuk{h), 

j,k=l 

(3.54) 

r22=n-3/2 ^ U,k{tl)Ujk{t2)Ujk{t3). 
j,k=l 

It follows from the Schwarz inequality and (2.16) that 

n / n \'^/'^/n \ 1/2 

hence, IT22I <n~-^/^. Besides, writing 

Tsi = n-^/^{U{ti)V{t2),V{t:i)), V{t) = n-^/\Uu{t), • . • , Unn{t)f, 

where by (2.16) ||y(t)|| < 1, \\U{t)\\ = 1, we conclude that |T2i| < n"i/^ 
hence, 

(3.55) IT2I <C2(t,x)/ni/2. 

This together with (3.40), (3.47), (3.51) and (3.52) imply that the r.h.s. of 
(3.43) is 0(n~^/^) as n ^ 00. We obtain (3.41), hence the assertion of the 
theorem. □ 

In fact, we have more (see [18] for a particular case of traces of resolvent). 
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Theorem 3.4. Theorem 3.3 remains valid if its condition (3.37) is re- 
placed by the Lindeberg type condition for the fourth moments of entries of 
W 

(3.56) lim L(^)(t) = Vr > 0, 

n— >oo 

where [cf (3.32)] 

1 " r 

(3.57) 4')(r) = 4 E / W'F^lXdW). 

Proof. Consider again the truncated matrix M'^ of (3.27). Since 

P{w ^w^}<j2 nw^k / 

j,k=l 

(3.58) 



= E / Fj-\dW)<r~'L(^Hr), 
we have, in view of (3.56), 

(3.59) hm E{e"'-^" - e^^^"- } = 0, 

where Nnri'^] = (p{M'^). Now it suffices to prove that if A/'n[(/'] = Tr(/3(M) 
is the hnear eigenvalue statistics of the GOE matrix Af , then [cf. (3.41)] 

(3.60) Rl{x) = E{e^^-^-[^l} - E{e^^-^"°['^]} 

vanishes after the hmit (3.34). To this end, we use again the interpolation 
matrix M'^(s) of (3.35), and get analogs of (3.43)-(3.46): 

K{x) = -^J^ds J t^{t)[Anr-Bnr]dt, 

(3.61) = ^ E{wj^^^^nr}, 



where now 



^nT = e°n^{s,x)UJk{s,t), 

(3.62) 

U^{s,t) = exp{itM^{s)}, enris,x) = exp{ixTTip{M^{s))°}. 
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Applying (3.6) with p = 3, we obtain [cf. (3.47)-(3.49)] 

3 

(3.63) Anr = Y.Tir + e3r. 

1=0 

(/-l)/2 n 

(3.64) Ti, = j^-^^^ Yl <+ijfcE{Z^jfc(5)^>nr}, ^ = 0,1,2,3, 

j,k=l 

and 

l^3r|<^ E sup |I),\(s)$„.|m(.)=m|- 

Since E{|W^J^^^|^} < Ty/n^i, then we have, in view of (3.50), 

(3.65) \e'ir\<Ci{t,x)T. 

Besides, it follows from (3.50) and (3.30) that we can replace Ti^ of (3.64) 
by Ti of (3.48) with of (3.62): 

(3.66) Tir = Ti + ri, 
where the error term r/ satisfies 

(/— 1)/2 

In I < /^(;+i)/2 E \^l+l,jk-^l+l,jkmD'jkis)'^nr}\ 
j,k=l 

(3.67) 

We have by (3.2) Tq = 0, Ti = Bn, and T2 and T3 satisfy (3.55) and (3.52), 
respectively. This together with (3.40), (3.43) and (3.65) imply (3.60) and 
complete the proof of the theorem. □ 



3.4. Central limit theorem in general case. Here we prove the CLT for 
linear eigenvalue statistics of the Wigner random matrix not assuming that 
the fourth cumulant of its entries is zero [see (3.39)]. We use the scheme 
of the proof of Theorem 2.2, based on the Gaussian differentiation formula 
(2.20) and the Poincare type "a priory" bound (2.26) for the variance of 
statistics. We have the extension of (2.20), given by (3.6). As for an analog 
of (2.26), it is given by the theorem below [see also (3.90)]. 

Theorem 3.5. Let M = n~'^/'^W he the Wigner matrix (3.1)-(3.3) sat- 
isfying (3.38) and (3.56), M'^ he corresponding truncated matrix (3.27), and 



(3.68) Unr{t) = TiU^{t), U^{t)=exp{itM^). 
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Then for any fixed r > 0, 

(3.69) Ysir{unr{t)} < a(M4)(l + |^|^)^ 

(3.70) Var{AA„, < a(/"4) (1 (1 + \t\^)\m\ dt 
where Cr{^J-4) depends only on fi4 and r. 

Proof. Note first that by the Schwarz inequahty for expectations and 
(2.60) we have 

Yar{J\fnr[v>]} = I I ^(tl)<^(t2)E{<^(tiX^(t2)} dtl dt2 

(3-71) 

and it suffices to get bound (3.69) for 

K = Var{n„,(t)}. 
Denoting Un{t) = TrexpjitM}, where M is the GOE matrix, we write 
Vn = E{n„(t)n°,(-t)} + E{(n„,(t) - 2„(t))2° (-t)} 

(3.72) 

+ ^{{unr{t) - 2„(t))«,(-t) - n° (-t))} = Ki+K2 + Ks. 
We have, by the Schwarz inequality and (2.65), 

(3.73) \Ki\ < V2w\t\V^^^, \K2\ < V2w\t\V^/'^ + 2w'^t^. 
To estimate K-^, we use the interpolating matrix (3.35) to write 

(3.74) Ks='jj\K^-B'^]ds, 
where 

with 

(3.76) = UJ,is,t)i<ri-t) - Ki-t)) 

and U'^{s,t) being defined in (3.62). 
Applying (3.6) with p = 2 and 

•^(W) := '^Wn\MJ^(s)={s/n)-^/^W+a-s)y2M.J 
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to every term 

E{t^t^^<,} = / HW)WF\f^{dW) 
of the sum in j4'„ (3.75), we obtain 

2 

(3.77) < = ^I); + e2r, 

1=0 

where T/^ is defined by (3.64) with of (3.76) instead of $„(s) of (3.62), 
and 

(3.78) 1^2.1 sup \nDU^mw)}\. 

Since 

E{I)jfc(sX} 

(3.79) =^{{ul,{-t)-ul{-t))D]k{s)U]„[s,t)} 

+ E (^) nD%{s){Unr{-t)-Un{-t))D\-\s)U],{s,t)^ 

and by (2.28), 

d 

(3.80) __Unr{t) = 2i(3,ktUJkit), 

jk 

the Schwarz inequality and (2.19) yield 

(3.81) \E{D^j,{s)K}\<Q{t){V^/' + l). 

Here and below we denote by Ci{t) an n-independent polynomial in \t\ of 
degree I. This and (3.30) imply [cf. (3.66) and (3.67)] 

(3.82) Ti = Tl + r[, Z = 0,l,2, 

where T/ is defined by (3.48) with of (3.76) instead of ^>„ of (3.46), 

= 0, K2jk = (1 + '^jfc)''^^, K3jfc = /is, and 

|r;|<.('-i)/2Q(t)r'-^L(f)(r)(yV2 + l). 
Taking in account (3.56), we have for sufficiently large n 

(3.83) \r[\ < s('-i)/2Q(t)^'-3(yV2 + i)_ 
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We see that = 0, and by applying (2.20) to B'^ of (3.75), we have T[ = B'^. 
Besides, since by (3.79) 



2 



j,k=i 



q=l ^ ' ^ j,k=l 

xD%''{s)UJ,{s,t)'jy 

then using the Schwarz inequahty and (2.19) to estimate the first term, and 
(3.80) and the argument leading to (3.55) to estimate the second term, we 
obtain 

\n\<c2mv^/'+i). 

It follows from the above for the integrand in (3.74) 

(3.84) K - B'J < \e2r\ + C2(t)(Ki/2 + i), 
where £2t is defined in (3.78) and we have, in view of (3.79), 

(3.85) \e2r\<^ E 5.^ + 14.1, 

j,k=l 



where 



S,k= sup |E{«,(-t)-n°(-t)) 

\W\<T./fi 



X -^jfc^/fc(^'*)lA./j^{s)={<i/n)i/2M/+(l-s)i/2Mj,,J"l 

and by (2.16), (2.19) and (3.80), 

(3.86) le'aJ < Csit) 

with C3(t) of (3.81). 

To estimate Sjk, we repeat again the above interpolating procedure, and 
obtain for every fixed pair {j, k} 



Sjk = Tr sup 



2 \W\<Ty^ 



»1 n 

/ ds,— y: E{(sr^/X^^^-(i-^i)"'^X.) 



^ *^nlMJJs)=(s/n)l/2vy+{l-s)l/2Mj.^} 
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where 

(3.87) K = U;g{si,t)D%UJk{s,t), \K\<C3{t). 

By the condition \W\ < T^/n and (3.87), two terms of the sum corresponding 
to Wpq = Wjk = W are bounded by C3(t) for every fixed r > 0. Hence, 
applying (3.6) and (2.20) to the rest of the terms, and using the notation 
J2'p^q for the sum with {p, q} / {j, k} and {p, q} ^ {A;, j}, we obtain 



5,fc<C4(t) + ^ sup f\A';,-B';\dsi 



(3.88) 
with 

1 / ^ 

and 

^ri = + '^P9)E{^P9('5l)^nlM'-Js)={s/n)l/2H/+(l-s)i/2M^-J' 

p,q 

where [cf. (3.64)] 

S^a-l)/2 , ^ ^ 
-^'^ ^ 71 (i+l)/2^ '^^+l,Pg^"L-^Pg(^l)*^nlM'-^(s)=(s/n)i/2H/+(l-s)i/2i/jfcJ' 

and 



7fe(«)=(V")''''l^+(l-s)i''2M^.fc,Af-(si)=Ml- 

Since |L>j,g(s)^>'^| < Q+3(t), then je'sVI < Ce{t). Besides, in view of (3.30), we 
have an analog of (3.82) and (3.83): 

r/; = T/' + rr, Z = 0,l,2, 

where an argument, similar to that leading to (3.55), implies 

\T^\<C5{t)n-'/^. 

We conclude that, for every r > 0, 

sup ['\A':,-B';^\dsi<C,it). 

Plugging this estimate in (3.88) and then in (3.85), we obtain in view of 
(3.86) that 

\e2r\ <C7{t). 
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This, (3.74) and (3.84) imply inequality \K^\ < C3{t)Vn^'^ + Csit), which 
together with (3.72) and (3.73) allow us to write the quadratic inequality 
for Vn: 

Vn<C3{t)V^/^ + Csit) 

valid for every fixed r > and any real t and implying (3.69). □ 

Remark 3.3. A similar but much simpler argument allows us to prove 
that if 

(3.89) 'u;6:=sup max E{\wi^^f}<oo, 

n l<j<k<n J 

then we have the bounds 

(3.90) Var{AA„ < C{we)(^J (1 + 1*1=^)1^(01 d?j^ 
and 

(3.91) VarK(i)} < C{wG)il + \tff, 

where C{wq) depends only on wq. The proof is based on the representation 

Va.r{un{t)} = B{un{t)ul{-t)} + E{(u„(t) - n„(t))n° (-t)}, 

the interpolation procedure, and the differentiation formula (3.6) with p = 4 
in the second term. 

Now we can prove the corresponding CLT. 

Theorem 3.6. Let M = n~^^'^W be the real symmetric Wigner matrix 
(3.1)-(3.3), satisfying (3.38) and (3.56), and (^crM^M be a test function 
whose Fourier transform (p satisfies (3.40). Then the centered linear eigen- 
value statistic [see (2.47)] converges in distribution to the Gaussian 
random variable of zero mean and variance 

(3.92) nv,bl=FooEM + ^(£%(M)^^g^<i.)^ 

where Vgoe[^] is given by (2.48), and K4 = — 3w'^ is the fourth cumulant 
of the off-diagonal entries ofW. 

Proof. Following the scheme of the proof of Theorems 2.2, we show 
that the limit Z{x) of characteristic functions Zn{x) = E{exp(ixA/'° [(/?])} 
satisfies (2.55) with Vgqe of (2.48) replaced by Vwig of (3.92). In view of 
(3.59), it suffices to find the limit as ?i ^ 00 of the characteristic functions 

(3.93) Z„^(x) =E{e„^(x)}, Cnrix) = ex.-p{ixN°^[ip\} 
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of the centered statistics M^ri'^] of truncated matrix M'^ = n~^^'^W'^ of 
(3.27), and then pass to the hmit r ^ 0. 

It is easy to see that formulas (2.57)-(2.64) with Unr and y„T-(x,t) = 
^{unrit) X e°^(x)} instead of n„ and Yn{x,t) of (2.61) and (2.63) are vahd 
in the Wigner case as well, and that (3.69) and (3.70) imply the analogs of 
(2.67) and (2.69) for Ynr- 

(3.94) \Ynr{x,t)\<Cl/^{fXi){l + \t\)^ 

and 

d 



(3.95) 



Ox 



<a(/U4)(y'(i+iti')i^(t)Miy 



where Cri^i) depends only on r and 

To prove the uniform boundedness of dYn{x,t)/dt [an analog of (2.68)], 
we note first that, by (3.1) and (3.68), 

o • ^ 

(3.96) j-Y^rix,t) = E{<,(t)e°,(x)} = ^ E ^{W^k'^^'^n}, 
where 

(3.97) $„ = C/;fc(t)e°,(x), \D'jk^n\<Ci{t,x), 0</<5 

[see (3.50)]. Treating the r.h.s. of (3.96) as A^r of (3.63) and applying (3.6) 
with p = 2, we obtain 

(3.98) -Yr,r{x,t) = '-^ {l + Sjk)B{D,k<^n}+0{l), 

j,k=l 

where the error term is bounded by C3(t, x) as n ^ 00 in view of (3.7), 
(3.55), (3.67) and (3.97). By using (2.17) and (2.72), we obtain for the first 
term of the r.h.s. of (3.98) 

itw'^n~^YnT{x,t) + iw'^ / ^In^^Unrit — ti)}YnT{x,ti) dti 

Jo 

+ iv? f ^{rT^UnAt - iiX,(ii)e°,(x)} dt^ 
Jo 

— 2w'^x J ti^(ii)E{n^^ti„T-(t + ti)e„T-(x)} (iti, 

where the first two terms are bounded in view of (3.94), the last term is 
bounded by 2ii;^|rc| / |ti||(^(ti)| dti, and the third term satisfies 

(3.99) \B{n-\nAt - ti)<,(ti)e°,(x)}| < 2Cr{^li)^/Hl + |ti|)^ 
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in view of (3.69). It follows then from (3.96)-(3.99) that, for any fixed r > 0, 



(3.100) 



r^^^r (^7 x] 



<C5{t,x). 



Thus, we have analogs of (2.67)~(2.69), implying that the sequence {l^nr} is 
bounded and equicontinues on any bounded set of M?. We will prove now 
that any uniformly convergent subsequence {1^;^} has the same limit 1^- 
To derive the limiting equation for Yr, we treat Ynr as Yn of (2.63), and 
applying first the Duhamel formula (2.14), write 

rt " 

Ynr{x,t) = ^ ^ E{T^^)^$„}dtl 
V ™ -'0 j ,,^-^ 

with of (3.97) [cf. (2.70)]. Then an argument, similar to that leading to 
(3.63)-(3.67) and based on (3.6) with p = 3, yields 

(3.101) Ynr{x, t)=if^ (jy.Ti + n) + e3r,n^ dh, 

where [cf. (3.64)] 

1 " 

(3.102) Ti= Yl ^i+ijk^{D'jk^n}, 1 = 0,1,2,3, 

^•^ j,k=i 

(3.103) \e3r,n\<rC4{t,x) 

and ri satisfies (3.67) with s = 1. Besides, Tq = 0, T2 satisfies (3.55), and the 
contribution to T3 due to the term 9w^6jk of K^jk = K4 — 9w^6jk [see (3.2) 
and (3.4)] is bounded by C{t,x)n~^ . This allows us to write an analog of 
(2.71) with additional term proportional to ^4: 

(3.104) y„,(:c,t) = r;2,„ + r;^„+f3r,n(x,t) + o(l), n^oo, 
where 

f* 1 " 

(3.105) T^2 n = iw^ - E (1 + 83k)^{Djk^n} dh, 

(3.106) T:^,n = i^4 / ^ E nD%^n}dh, 

•^0 j,k=l 

(3.107) 83r,n{t,x)= [ e3r,n{tl,x)dti 

Jo 

and for any r > the reminder term o(l) in (3.104) vanishes as n — > 00 
uniformly on any compact set of {i > 0,a; S ffi}. The term T^2 „ of (3.105) 
has the same form as the r.h.s. of (2.71) of the GOE case. Since the argument. 
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leading from (2.71) to (2.76)-(2.78) does not use the Gaussian form of Wjk 
in (2.71), it is applicable in our case as well and yields 

'^w^ n = -2^«^ / dtl [ VnT{t2)YnT{x,tl - 12) dt2 

Jo Jo 

+ xZm-{x)Am-(t) — r^T-(x, t), 

where Vnr = n^^^iunr}, and Anr and rnr are given by (2.77) and (2.78) 
with the GOE matrix M replaced by the truncated Wigner matrix M'^ . 
Now it follows from the Schwarz inequality, (3.69) and (3.40) that, for any 
T > 0, the reminder rnr{x, t) vanishes as n ^ oo uniformly on any compact of 
{t<0,x£ M}. Besides, in view of |U„t-| < 1 and (3.58), linin^ooivnT — Vn) = 0, 
Vr > 0, and Theorem 2.1 yields that, for any r > 0, the sequences {vut} and 
{74„T-} converge uniformly as n — > oo on any finite interval of R to v{t) and 
A{t) of (2.80) and (2.81). It follows also from (3.93), Theorem 3.5 and (3.40) 
that 

\ZUx)\ <\x\VSir^/^Mnrm 



< \x\Cr{H4) /(I + l*l'^)l<?(OI dt < OO. 



Hence, the sequence {.^nr}n>o is compact for any r > 0. Denoting the con- 
tinuous limit of some its subsequence {ZniT}n>o by Zr, we have for any 
r > uniformly on any compact set of {t > 0, x G M} 

lim T^2^ =-2w^ f dtl r v{t2)Yr{x,ti-t2)dt2 + xZr{x)A{t). 

ni^oo > ' Jo Jo 

(3.108) 

Consider now the term T^^ ^ of (3.106) and note first that, in view of (3.58) 
and (3.97), we can replace T^^ „ by 

(3.109) T,,,n = ^ E nD^k{Ujk{ti)eUxm dt 



1 

j,k=i 

4r{4)/ 



with the error bounded by C4^{t,x)T Ln (t). It follows now from (2.17), 
(2.72), (3.54) and (3.55) that the contribution to T^^^n due to any term of 



J2 D%iU,kit)el{x)), 
j,k=i 

containing at least one off-diagonal element Ujk, is bounded by C3{t,x)n~^ . 
Thus, we are left with terms, containing only diagonal elements of U. These 
terms arise from e°^{x)D%U jk{t) and 3Dji^Ujk{t)D'^f,e'^{x) in the above sum, 



and, by (2.17) and (2.72), their contributions to T^^^^n are 

n 

(3.110) ^ Yl / E{(f/jj * Ujj * Ukk * Ukk){ti)el{x)} dtl 
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and 

(3.111) '^J2 dtJ t2^{t2)B{{U,j *Ukk){tl){Uj, *Ukk){t2)en{x)} dt2, 

" j,k=l 

where we omitted because the corresponding error term is 0(n~^). It is 
easy to see that the entries of U appear in (3.110) and (3.111) in the form 

(3.112) B{vn{tut2)vn{t3,U)el{x)} 
and 

(3.113) ^{Vn{tl,t2)Vn{t3,h)enix)}, 

where 

n 

(3.114) Vnitl,t2) =n~^J2^j,ih)U,jit2). 

i=i 

Since \Ujj{t)\ < l,t G M, we have 

(3.115) |^n(tl,t2)| <1. 

This, the inequaUty |e° (x)| < 2, and the general inequality 

(3.116) E{|(66)°|} < 2cE{|^i°|} + 2cE{|^2°l}, 

where 2 = Ci,2 — E{^i^2}) and ^1^2 sue random variables such that |^i,2| ^ c 
allow us to write for (3.112) 

|EK(ti,t2)t^n(t3,i4)e°(x)}| <4E{|t;°(ti,t2)|} + 4E{|<(t3,t,4)|}. 
By Lemma 3.1 below, we have 

(3.117) E{|<(ti,t2)|} < C3ih,t2)n-y^. 

Thus, (3.112) vanishes as n — > 00 uniformly in t and x, varying in any com- 
pact set of . 

Expression (3.113) can be written as the sum of (3.112) and 

(3.118) B{Vn{h,t2)Vn{t3:U)}E{en{x)}=B{Vn{tut2)Vn{t3:U)}Zn{x). 

It is follows from (3.115) and (3.117) that E{vn{ti,t2)vn{t3,ti)} can be 
written as the product Vn{ti,t2)vn{t3,t4) up to an error term bounded by 
C3(ti,t2)^/2C3(t3,t4)^/'n-i/4, where 

(3.119) Vn{ti,t2)=B{Vn{tl,t2)}. 

In addition, we have, by Lemma 3.1 below, 

(3.120) vnituh) = v{ti)v{t2) + 0(1), 

where v is given by (2.80) and o(l) is bounded by (73(^1, ^2)"-"^''^. 
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We conclude from the above that the contribution of (3.110) to TK4,n; 
vanishes as n/ — > oo uniformly in t and x, varying in any compact set of 
{t>0,x£ M}, while in (3.111) we can replace Ujj and Ukk by v. As a result, 
we obtain 



(3.121) 



T^^ := lim T^,^ „ 



ixZT-{x)K4 / {v*v){ti)dti t2(p{t2){v *v){t2)dt2, 







uniformly on any compact of {t>0,xGM}. 

In view of (2.83) and Proposition 2.1(iii), we have 

{v*vm = -^l^e^'^f\z)dz. 

The integral over L can be replaced by that over the cut [—2w, 2w] of 
\J — Aw"^ in (2.36) and we obtain that 

(3.122) („*^)(t) = __!_ l'\'^y^4w^-^,^dt, 
or, integrating by parts, 

1 i-2w 9, ,,2 _ .,2 

(3.123) e-^^^^d,. 



TTtW^ J-2w y^4w^ - fJ? 

Now the Parseval equation implies that 

(3.124) / tif{t){v * v){t) dt = ^ ififi) ^ d^i =■ B, 
thus, 

T^^ = iBI{t)xZr{x), 

where 

(3.125) I{t)= f {v*v){ti)dti. 

Jo 

Besides, it follows from (3.104) and the convergence of sequences {Yfur}, 
{r^2 „J and {rj_j,„J [(3.108) and (3.121)] that the limit 

(3.126) S3rit,x) = lim S3r,n,{t,x) 

rii — ^oo 

exists uniformly on any compact of {i > 0,x G M}, and we have by (3.103) 
and (3.107) 

(3.127) \£3r{t,x)\<TC5it,x). 
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This and (3.108) allow us to pass to the hmit n; — > oo in (3.104) and to 
obtain the integral equation 

Yr{x,t)+2w^ f dti f ' v{ti-t2)Yr{x,t2)dt2 



(3.128) "° 

= xZr{x)[A{t) + iKiBI{t)] +£3r{t,x). 

The l.h.s. of (3.128) coincides with that of (2.82) and the r.h.s. of (3.128) is 
equal to that of (2.82) plus two more terms. Thus, the solution of (3.128) 
is equal to the r.h.s. of (2.86) plus two more terms, the contributions of the 
additional terms in the r.h.s. of (3.128). The r.h.s. of (2.86) leads to the 
first term in (3.92) [see (2.87) and the subsequent argument]. To find the 
contribution to (3.92) of the second term of the r.h.s. of (3.128), we use the 
r.h.s. of (2.11) with R{t) = iKixZ^{x)BI{t), Ti of (2.85), and (3.122). This 
leads to the term 

ixZr{x)B f^^ e»fA(2w^ - A^) 

in the solution of (3.128), where we used the relations 

dX = 0, 



2w \/4to2 - A2(A - fi) 



r2w V4u;2 - A2 ,^ , , ^ 

/ — ^ — flA = — vr/i, \fj.\<2w. 

J~2w (A - ii) 

Then the limiting form of (2.62) and (2.57) yield the expression k^^B"^ /2, 
that is, the second term of (3.92). 

Let us consider the contribution C3T-(t,x) of the third term of the r.h.s. 
of (3.128), which is given by the r.h.s. of (2.11) with R{t) =£3r{i,x) and Ti 
of (2.85). Integrating by parts, we obtain 

(3.129) C3r{t,x)=Ti{0)£3r{t,x)+ f T[{t-t{)83r{tl,x)dh. 

Jo 

We have also 

Ti(i) = -Jo{2wt), T[{t) = 2wJi{2wt), 

where Jq and Ji are the corresponding Bessel functions, so that |Ti(t)| < 1, 
|r{(i)| <2w, and 

T[{t) = ^f^sm{2wt - 7r/4)(l + 0(t-3/2)), t ^ oo. 

V TTt 

By using this, (3.129) and (3.127), it can be shown that 

(3.130) \C3rit,x)\<TC5it,x). 
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Now, the limiting form of (2.62) implies 

(3.131) Z'^{x) = -xVw.^Zr{x) + Dsrix), 
where 

D3rix)=i J (p{t)C3rit,x)dt 

and, in view of (3.130) 

(3.132) |^3r(a;)| < tC4{x) /(I + \tf)\'f{t)\ dt, 



and C4 is n- and r-independent polynomial in \x\ of degree 4. 
Since Zt-(O) = 1, we can replace (3.131) by 

Jo 

and then (3.40) and (3.132) imply that 

limZ^(x)=e-^^'^^'/2 
hence the assertion of theorem. □ 

Remark 3.4. (1) The proof of the CLT under condition (3.89) is much 
simpler, because it does not use the truncation procedure. 

(2) Another expression for the limiting variance of linear eigenvalue statis- 
tics is obtained in [2]. In fact, the paper deals with the more general class of 
random matrices that the authors called the band matrices and that includes 
the sample covariance matrices with uncorrelated entries of data matrices 
of Section 4 below. Thus, a rather general formula for the variance of linear 
eigenvalue statistics obtained in [2] reduces to formulas (4.28) and (4.65) 
below. 

It remains to prove the following: 

Lemma 3.1. We have under the assumptions of Theorem 3.6 
(3.133) Yar{vn{ti,t2)] < C3(ti, ^2)/^^/^ 

and 

C3{ti,t2) 



'E,{Vn{tl,t2)} = Vn{tl)Vn{t2)+rn{tl,t2), |r'„(tl,t2)| < 

(3.134) 

where Cp is a polynomial in \tj\, j = 1,2, of degree p with positive coefficients. 
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Proof. We denote again M = n'^/^W the GOE matrix (2.23), U{t) 
e^**^ Vniti,t2) = Ej=i Ujj{ti)Ujj{t2), and write 

Var{vn{ti,t2)} = E{Vn{tl,t2)v°^{-tl,-t2)} 

(3.135) +E{{Vn(.tl,t2)-Vnitl,t2))v°^{-h,-t2)} 

= : i?i + i?2- 

The Poincare inequahty (2.17), (2.21) and (2.24) ahow to obtain 

(3.136) Var{u„(ti,t2)} <4^i;2(^? + ^i)n-2 
and 

(3.137) VariUjjit)} < 2w'^t'^n~\ 
hence, 

(3.138) \Ri\<Ci{ti,t2)n-^. 

To estimate R2, we use again the interpolation matrix (3.20) and write 
R2=l —-Y,^{Un{s,ti)Ujj{s,t2K{-ti,-t2)}ds 

■ I E ms-'^'w^^-{i-s)-'^'w,k)^{s,h,t2)}ds, 



pi n 

2n3/2 



where 

^'(S, ti, ta) = ^^n(-il> -i2)(tlf/j/fc(s, tl)f/ji(s, t2) + t2f/jfc(s, t2)Ujj{s, ti)) 

and f/(s,t) is defined in (3.46). We have by (2.20), (3.6) with p = 1, and 
(2.19) 



\R2\<C^{tiM)n 



-1/2 



[cf. (3.24)]. This, (3.138) and (3.135) yield (3.133). 
To prove (3.134), we write 

(3.139) ¥.{Vn{tl,t2)] = E{Vn{tl,t2)] + E{Vn{tl,t2) -Vn{tl,t2)], 

where the second term similar to term R2 above is modulo bounded by C3 x 
(ti,t2)n-V2. 

It follows from the orthogonal invariance of the GOE that E{i7jj(t)} = 
U„(t), so that 

(3.140) E{Vn{tl,t2)] =^n,(tl)^n(t2) +E|n-lX^^C/,-,(ti)C/°^.(i2)|, 
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where by (3.137) the second term is modulo bounded by Ci{ti,t2)n ^. Be- 
sides, by using interpolation procedure, we can show that 

(3.141) mvnit) - Vnim < C3{t)n-^/\ 

This, (3.139) and (3.140) yield (3.134). □ 

4. Sample covariance matrices. 

4.1. Generalities. We again confine ourselves to the real symmetric ma- 
trices. Thus, we consider in this section n x n real symmetric matrices 

(4.1) M = y^y, Y = n~^/'^X, 

where X = {X^'^^}^'^^^ is the m x n real random matrix with the distri- 
bution 

m n 

(4.2) PmnidX) = n n Fi7'''\dX^j)^ 

a=lj=l 

satisfying 

I XFi7'")(dX) =0, I X^Fi^'''\dX) = a\ 
In other words, the entries {Afj™'"'''}"^^-|^ of Af of (4.1) are 

m 
a=l 

where X^^'"^ G R, a = 1, . . . , m, j = 1, . . . , n, are independent random vari- 
ables such that 

(4.3) E{47''^)} = 0, E{47''^)4'^'")} = ^'^P^^^^''- 
A particular case of (4.1)-(4.3), 

(4.4) M = y^y, y = n"^/2^, 

where the entries of X = {Xq.jI^j'Li ^-re i.i.d. Gaussian random variables 
satisfying (4.3), that is, 

m n 

(4.5) P(dX) = Z-^i exp{- TTX^X/2a^] J] 11 ^^-i' 

Q = l j = l 

is closely related to the null (white) case of the Wishart random matrix of 
statistics (see [21], Section 3.2). The difference is in the factor m~'^/'^ instead 
of n^^/^ in (4.4). In what follows, to simplify the notation, we will often omit 



46 



A. LYTOVA AND L. PASTUR 



the superscript {m,n), and the sums over the Latin indexes wih be from 1 
to n, and the sums over the Greek indexes will be from 1 to m. 

We present first an analog of the law of large numbers for the sample 
covariance matrices. 



Theorem 4.1. Let M be the real symmetric sample covariance matrix 
(4- 1)-(4-3)- Assume that for any r > 

,4.) j^j:/ x=^r'(«)-o 



as 



(4.7) n^oo, m^oo, m/n ^ c € [0, oo), 

and that {Xaj}^'j^=i are defined on the same probability space for all m,n £ 
N. Then for any bounded continuous c^iM— >C, we have with probability 1 

(4.8) lim ^ n"W„M = / ^{X)NMp{d\), 
where Nn['-p\ is defined in (1-1), and 



(4.9) NMp{d\) = (1 - c)+5o(A) d\ + (2^a^\)~^ ^ {{\- a^){a+ - X))^d\ 
with a± = {1 ^ y/c)'^ , and 3;+ = max(x,0). 

We refer the reader to [4] and [12] for results and references concerning 
this assertion that dates back to [20]. Here we outline a weaker version of 
the theorem on the convergence in mean in (4.8), basing it on the same ideas 
as in Theorems 2.1-3.2. We will need this assertion as well as the method 
of its proof. 

We start from the Gaussian case, that is, the Wishart random matrices 
(4.4) and (4.5), and follow essentially the proof of Theorem 2.1. Introduce 
[cf. (2.37)] 

(4.10) /„(z) = E{ff„(z)}, gn{z)=n-^TrG{z), G{z) = {M ■ 



z)-\ 



By using again the resolvent identity (2.40) and the Gaussian differenti- 
ation formula (2.20), we obtain for /„ of (2.37) 

fn{z) = -i + -i_ ^E{X„fc(yG)„fc(z)} 

a,k 

1 a2 



z zn 
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We have from (2.40) and (4.4) [cf. (2.42)] 

(4.11) DakGjk = -{YG)akGjk — {yG)ajGkk, 

hence, by (2.44) [cf. (2.43)], 

fn{z) = -- + ^fn{z) - -E{gn{z)n-' TrG{z)M} 

(4.12) 

-—E{n-'^TrG^{z)M} 
z 

or, after using the identity 

(4.13) G{z)M = zG{z) + I, 

we get 

1 

fn{z) = + —Cnfn{z) - — E{5f„ (z) (z5„ (z) + 1)} 

a2 



— E{n-iTVG(z)(zG(z) + l)} 
nz 

with 

(4.14) Cn = nn/n. 

We need now the Poincare type inequahties for the Wishart matrices (4.4) 
and (4.5): 



(4.15) Var{AA„[(^]} < Aa^E{n-^ Ti ip' {M)ip' {M)M} 

(4.16) <4a'^c„sup|v3'(A)|2 



[cf. (2.25) and (2.26)]. They can be easily derived from (2.21) by using the 
formulas E{n~^ TrM} = a^c„ [see (4.3)] and 

(4.17) D^kUji{t) = mYU)o.j * Uki){t) + {{YUU * C/,fc)(t)), 

(4.18) D^kTrip{M) = 2{Y^'{M)U 

[cf. (2.17) and (2.28)]. By applying (4.16) with 93(A) =n~^{X - z)-\ we 
obtain [cf. (2.45)] 

4(7^r 

(4-19) Var{5„,(z)} < 

Thus, (2.40) and (4.12) allow us to write 

fn{z) = -z~^ - a^z~^{l - Cn)fn{z) - fl{z) + 
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where 

, , Aa^Cn 1c? 

\rn\ < ni^. I. + 



nlQz? 



This and the hmit (4.7) yield an analog of (2.46): 

(4.20) za''fMp{z) + iz + a\l-c))fMp{z) + 1 = 0, Qz^O, 
and since QfMp{z)Qz > 0, we obtain [cf. (2.36)] 

(4.21) fMp{z) = {^{z - a™,)2 - Aa^c - (z + 0^(1 - c)))/2a'z, 

where the branch of the square root is fixed by the asymptotic form z + 0(1), 
z — > oo, and 

(4.22) am = a^{c + l). 
This and inversion formula 

Nmp{^) = lim vr^i / 9/a/p(A + iO) d\ 

where the endpoints of A are not the atoms of Nmp, lead to (4.9). 
The next step is to prove an analog of Theorem 3.1, assuming that 

(4.23) 03:= sup max E{|xi7'"^|^} < 00. 

To this end, we use again an "interpolation" matrix [cf. (3.20)] 

(4.24) M{s)=Y^{s)Y{s), Y{s) = s'^/^Y + {1 - s)^/^Y, sG[0,1], 

where Y and Y are defined in (4.1)-(4.5). We have, with the same notation 
as in Theorem 3.1, 

fn{z) - fn{z) 

(4.25) = r -^Bhi-^TrG(s,z)}ds 

Jo OS 

-3/2 /'^E{(s-i/2x„fc - (1 - sr'?X^,){Y{s)G'U}ds. 
Jo ^ 



-n 



1 

1 m,n 



Since {Xaj}^'j^i are independent random variables satisfying (4.3) and 
(4.23), and {X^.j}'^'] are i.i.d. Gaussian random variables also satisfying 
(4.3), we apply the general differentiation formula (3.6) with <1> = {Y{s)G')oLj 
and p = 1 to the contribution of the first term in the parentheses of (4.25) 
and the Gaussian differentiation formula to the contribution of the second 
term. As it was already several times in the case of the Wigner matrices 
(see, e.g.. Corollary 3.1 and Theorem 3.1), the term with the first derivative 
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of the general differentiation formula is canceled by the expression resulting 
from the Gaussian differentiation formula, and we are left with [cf. (3.24)] 

fn{z) - fn{z) = [ V^ei(s)ds, 

Jo 

where 

|ei(s)|<4^5] sup \Dlj{YG'U\, D^k = d/dYak, 

where G = {Y'^Y - zy\ and M m,n is the set of m x n real matrices. It 
suffices to find an 0(1) bound for D'^^{YG')a,k- Since (YG)cik is analytic in 
z, '^z 7^ 0, then the bound for (YG')ak follows from that for {YG)ak and 
the Cauchy bound for derivatives of analytic function. By using (4.11) and 
a little algebra, we obtain 

Dl^iYGU = -QGkk{YGU + QGkk{yG)^k{YGY^)^^ + 2{YG)l^. 

It follows from (2.40) that \Gkk\ < Next, if G = {Y^Y - zY"" and^G = 

(yy^ - z)-i, then YG = GY , and {YGY'^)^,^ = {GYY'^)^,^ = (1 + zG)aa 
[see (4.13)], thus, 

(4.26) \{YGY^)^a\<l + \z\\'^z\-^. 
Furthermore, it follows from the Schwarz inequality that 

(4.27) |(yG)„fc| < (G*y^yG)^// < ((l + \z\\'^z\-^)/\'^z\)^'\ 

Thus, D^i^(YG)ak is bounded uniformly in 1 < a < m, 1 < k <n, all m and 
n, and z, varying in a compact set K C C \ M, and 

|ei(s)| <G_ft:n~"^/^, n ^ oo,m — > c>o,m/n — > c G [0,oo), 

where Gk < oo depends only on K C C \ M. 

In fact, a bit more tedious algebra and (4.26) and (4.27) show that, for 
every 1 < a <m, 1 <k <n {YG)ak{z) is real analytic in every y^j, 1 < (3 < 
m, I < j <n and Qz 7^ 0. Hence, all derivatives dl^j^{YG)ctk{z), / = 0, 1, . . . , 
are bounded by Gik, zEK cC\R [cf. (3.25)]. 

This proves (4.9) under condition (4.23), that is, an analog of Theorem 3.1. 
To prove (4.9) under condition (4.6), we have to use the truncation procedure 
analogous to that of the proof of Theorem (3.2) and bounds (2.40), (4.26) 
and (4.27). 
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4.2. Central limit theorem for linear eigenvalue statistics of the Wishart 
ensemble. Following our scheme of the presentation in the case of the 
Wigner matrices, we start from the central limit theorem for the sample 
covariance matrices with Gaussian entries, that is, from the Wishart ensem- 
ble (4.4) and (4.5). We confine ourselves to the case c > 1. 

Theorem 4.2. Let J\fn[f] be a linear eigenvalue statistic of the Wishart 
matrix (4-4) (4-5)j corresponding to a bounded function (/?:M— >]R with 
bounded derivative. Then the centered random variable M°[(p] [see (2.47)] 
converges in distribution asn— >oo,m— >oo,m/n— >c>l to the Gaussian 
random variable of zero mean and variance 



(4.28) X ^ 4a^c-(Ai-a^)(A2-a^) 

^ ^ V4a4c - (Ai - a„) V4a4c - (A2 - a^r 

where Aip is defined in (2.49), a± = a^(lib y/c)'^ and am is defined in (4.22). 



Proof. We follow the scheme of the proof of Theorem 2.2. Namely, 
assume first that admits the Fourier transform (p [see (2.53)], satisfying 
(2.54). We have, similarly to the proof of Theorem 2.2, the relations (2.57)- 
(2.64). It follows also from (4.16) with (/9(A) = e^*^ that [cf. (2.65)] 

(4.29) Var{u„(t)} < AaH'^Cn, 
thus [cf. (2.67)], 

(4.30) \Yn{x,t)\ = \B{ul{t)en{x)}\<\s.r^'\un{t)}<2a^\t\c]l\ 
Likewise, we have the bound 

(4.31) \dYn{x,t)/dx\ <4a^v^sup|(^'(A)|, 

following from (2.21) and (4.29) [cf. (2.69)], and the bound 

(4.32) \dYn{x,t)/dt\ < 20^^(1 + Ca^t2)i/2 

with C depending only on c„, following from (4.5) and (4.15) [cf. (2.68)]. 
Hence, the sequence {1^} is bounded and equicontinuous on any finite set of 
M^. We will prove now that any uniformly converging subsequence of {1^} 
has the same limit y, leading to (2.55), hence to (2.51) and (2.52) with 
Vwish[¥'] instead of Vgoe[v3]. Applying the Duhamel formula (2.14), (2.20), 
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(4.17) and (4.18), we obtain 

Yn{x,t) = ia^Cn [ Yn{x,ti)dti-a^n~^ [ dti / ' E{Tr M[/(ti)e° (x)} 
Jo Jo Jo 

-a^n~^ f dti /*'E{TrMC/(ti -t2)TrC/(t2)e°(x)}dt2 
Jo Jo 



or 



(4.33) 



10 JO 

la^xrT^ C E{Tr (/7'(M)MC/(ti)e„(j;)} dtx 



y„,(a;,t) = m^(c„ - 1) / Yn{x,t\)dtx 
Jo 

+ ia'^n-^ f B{u'^{ti)el{x)}tidti 
Jo 



+ ia^n ^ I 'Ei{un{t - ti)un{ti)en{x)} dti 
Jo 

+ 2ia^xn-^B{Tiip'{M){U{t) - l)e„(x)}, 
where we used the formulas TrMJ7(t) = —iu'^{t) and 

dti I ' E{<(ti - t2)Un{t2)el{x)] dt2 

Jo 

'B{{un{t - h) - n)un{ti)el{x)} dti. 
This and an analog of (2.74) and (2.75) yield an analog of (2.76), 

Yn{x,t) -ia^{cn-l) f Yn{x,ti)dti -2ia^ [ Vn{t - ti)Yn{x,ti) dti 



(4.34) 

= 2ia^xZnix) [ ip'{X){e''^-l)E{NnidX)} + rn{x,t), 



where now Vn = n ^E{n„} and 

rn{x,t) = ia^n~^ [ {Yn{x,t) -Yn{x,ti)) dti 



(4.35) +ia^n^^ T E{u° (t - tiX(ti)e° (x)} dii 

Jo 

- 2c? xn-^ J 9ip{9){Yn{x, t + 9)- Yn{x, 9)) dO. 

It follows from (2.54), (4.29) and (4.30) that rn{x,t) = 0{n~'^) uniformly 
in (x,t), varying in a compact set K C {x £M,t > 0}. This and Theorem 
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4.1 imply that the hmit of every uniformly converging subsequence of {Yn} 
solves the equation [of. (2.82)] 

Y{x,t)-ia^{c-l) f Y{x,ti)dti 

,2 







-2ia^ / VMp{t-ti)Y{x,ti)dti = xZ{x)A{t), 
Jo 

where [cf. (2.80)] 

(4.36) VMp{t) := lim U„(t) = — ^ / e^'^JAa^c - (A - a„,)2A^i dX, 

and [cf. (2.81)] 

A{t) = 2c?i [ (^'(A)(e**^ - l)NMp{dX) 



(4.37) 

1 /■* /■'^+ 



^ dh £^ e^*i^(^'(A)y^4a4c - (A - 0^)2 dX. 



Now an argument similar to that leading from (2.82) to (2.86) and based on 
Proposition 2.1 and the formula 

(4.38) VMP = Imp 

yields 

Y{x,t) = '-^^£^ if'{X)^Aa^c-{ix-amYdX 

(4.39) 



V4a4c-(/u-a™)2(^ - A) 



d^. 



Using this in (2.62), we obtain an analog of (2.87), and then an analog of 
(2.55) via (2.58) with Vwish of (4.28) instead of Vgoe, that is, an equation for 
the limiting characteristic function. Since the equation is uniquely soluble, 
we have finally 



Z{x) = e 



that is, the assertion of the theorem under condition (2.54). The general 
case of bounded test functions with bounded derivative can be obtained via 
an approximation procedure analogous to that of the end of the proof of 
Theorem 2.2 and based on (4.16). □ 

Remark 4.1. (1) The proof of Theorem 4.2 can be easily modified to 
obtain an analogous assertion for the Laguerre ensemble of Hermitian matri- 
ces M = n~^X* X, where the complex m x n matrix X has the probability 
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distribution [cf. (4.5)] 

m n 

P{dX) = Z„^2 exp{- TrX*X/a'^} Y[ Y[ d^Xo^j dQXo^j. 

0=1 j=i 

The result is given by Theorem 4.2, in which Vwish is replaced by VLag = 

Vwish/2. 

(2) It follows from the representation of the density p„ of E{A^„} via the 
Laguerre polynomials that (see [19], Chapters 6 and 7) 

for finite c and C and A sufficiently big. This bound and the approximation 
procedure of the end of proof of Theorem 2.2 allows us to extend the theorem 
to test functions whose derivative grows as Cie'^^'^ for some ci > and 
Ci < 00. 

4.3. Central limit theorem for linear eigenvalue statistics of sample co- 
variance matrices: the case of zero excess of entries. We prove here an 
analog of Theorem 3.4 for the sample covariance matrices. 

Theorem 4.3. Let M he the sample covariance matrix (^.l)-(4-3). As- 
sume the following: 

(i) the third and fourth moments of entries do not depend on j, k, m 
and n: 

(4.40) = E{ {xj-'-^f}, /.4 = E{ (Xi7'") )n; 

(ii) for any r > 0, 

(4.41) (^) ..= n-^Y. [ X^Fi-'-\dX)^0 

as oo, m— >-oo, m/n^ c£ [1, oo); 

(iii) the fourth cumulant of entries is zero: 

(4.42) Ki = ^ii- ^a'^ = 0. 

Let yjiM— >R he a test function whose Fourier transform satisfies (3.40). 

Then the corresponding centered linear eigenvalue statistic N^lip] con- 
verges in distrihution to the Gaussian random variable of zero mean and 
variance FwishM of (4.28). 

Proof. We follow the scheme of the proof of Theorem 3.4. Thus, in 
view of Theorem 4.2, it suffices to prove that if subsequently 

(4.43) m,n— >oo, m/re — > c G [1, oo) and r ^ 0, 
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then [cf. (3.60)] 

(4.44) R"mn{x) = Ele"^-!"^]} - E{e*^^"[^l} ^ 0, 

where Mnri^] is a hnear eigenvalue statistic corresponding to the truncated 
matrix [cf. (3.27)] 

(4.45) 

= {^l, = signXi-'") max{|47'")|, rnV2}}--'l^ 

and the statistic A/'n[97] corresponds to the Wishart matrix Y^Y of (4.4). By 
using interpolating matrix [cf. (3.35) and (4.24)] 

M^(s) = y^^(s)y^(s), 

(4.46) 

y-(s) = si/2y- + (i-s)V2y^ sG[0,l], 
we have [cf. (3.43)-(3.46)] 

(4.47) ^mn(^) = -^^^ *^(t)i^n " ^n] dt, 
where now 

with 

(4.48) ^ak{s) = e° (s, x){Y^{s)U{s, t)U 

and e„(s,x) and U{s,t) are defined in (3.62) in which M"^ is given by (4.45) 
and (4.46). We have, by (2.20), 

Bn = -Y.E{Dakis)^ak{s)}, Dakis) = d/dY;k{s) 
n , 

and, by (3.6) with p = 3 [cf. (3.63)-(3.65)], 

3 

An = '^Tir +e3T, 
1=0 

where now 

(l-l)/2 

(4-49) = E4+l,»fcE{-PL(^)^afc(^)}, ^ = 0,1, 2, 3, 

Q,fe 

ak cumulant of X'^^ , and 

C T 
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in view of E{|X^^,|^} < T^/nfli. 

In what follows we omit s and denote Dak = Dak{s), U{t) = U{t,s), etc. 
Let us prove the uniform boundedness of derivatives Eji^^^^afc}, / < 4, 
that will allow us to obtain analogs of (3.65)"(3.67). We have, analogously 
to (4.17), 

(4.50) DakUjkit) = im^UU * U,k){t) + {{YW)aj * UkkXt)), 

(4.51) Dakenix) = -2xen{x) J 9(p{9){Y^U)ak{0) d9, 

Dak{Y^UU{t) = Ukk{t) + i{{{Y^UY^^)aa * UkkM 

(4.52) 

+ {{Y^U)ak*{Y^U)akm), 

(4.53) Dak{Y^UY^^)aa{t) = 2{YW)ak{t) + 2i{{Y-UY-^)aa * {Y^UUm. 
Since by (2.16) 

(4.54) \{Y^UU{t)\ < (EiYajfY', \{YWY^^U{t)\<Y.iyaj?, 

j j 

then, iterating (4.50)-(4.53), we have 

/ N (« + l)/2 

(4.55) \DU<^>ak\<Q{t,x)r£{Y^jf] 

j 

and by (4.46), 



l>0. 



Now the Holder inequality implies the bound 

j j 
and analogous bounds for {Xaj}, thus, 
(4.56) \E{D'ak^ak}\<Ci{t), l<3. 

In the case where I = 4 a similar argument and (4.55) yield 



sup |E{D^fe$«fe| }| 



<C4(t,x), 



aj I 
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where we took into account that we have, by the Holder inequahty and 
condition \ < r^/n, 



(4.57) 



3\ 5/6 

aj 



j j 3 j 

< n-5/2(^2^2^^ ^ 3nV4a^ + n^a^^^ < C < oo 
with n-independent C. We conchide that 

(4.58) \e3r\<C4it,x)T 

[cf. (3.65)]. Besides, (4.56) and an analog of (3.30) for Xaj allow us to obtain 
for Tir of (4.49) an analog of (3.66) and (3.67): 

(4.59) Tir = Ti + ri, 
where now 

g{l-l)/2 

(4.60) Ti = J^^^(J^ E ^^i+i,ak^{DUs)^ak}, 

(4.61) \n\<Ci{t,xy-^L(^^{T). 

We have Tq = = (recall that Ki^ak = «^4,ofc = 0)) Ti = Bn, and, in view 
of Lemma 4.1 below, T2 = o(l) [cf. (3.55)]. Hence, 

An=Bn + e3r + 0{l), 

where the error term is a polynomial in |t| and \x\ of degree 3 at most 
that vanishes as m,n ^ 00, m/n ^ c uniformly in t and x varying in a 
compact set K c{t>0,x £ R}. This, (3.40), (4.47) and (4.58) imply (4.44) 
and complete the proof of the theorem. □ 

Remark 4.2. A similar argument leads to the proof of the CLT for lin- 
ear eigenvalue statistics of Hermitian analogs of (4.1)-(4.3), satisfying (1.5). 
The variance of the corresponding Gaussian law is V\vish/2, where Vwish is 
given by (4.28). For real analytic test functions this formula is a particular 
case of the variance, obtained in [5] for random matrices n~^X*TX, where 
X is a complex matrix with i.i.d. entries, satisfying (1.5) and (4.3), and T 
is a certain Hermitian matrix. 
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Lemma 4.1. We have under the conditions of Theorem 4-3 
(4.62) T2 = ^ E{Dlde°n{^){Y^U{t)U)} = o(l) 

as m,n^ oo, m/n c. 



a,k 



Proof. By using (4.50)-(4.53), it can be shown that the assertion wih 
fohow from 



r2p = o(l), p= 1,2,3, 



with 



72i=n-2^E{(X";7)«fc(fi)}, 



T22 = n-''Y.^{{X-U)^k{ti){X^UU{t2){X^U)^k{h)}, 

a,k 

T23 = n-''Y.^{{X^UX^^)^^{ti){X^U)ak{t2)} 

a,k 

[cf. (3.54)]. The Schwarz inequahty, (2.16) and (4.3) yield 
|r2i|=n-2|E{^(^C/,fc(tO) (E^«.)} 



j k 



and [see also (4.57)] 



h, 32,33 " 



<n-/2EV^{E(E(^.^.)0' + 3 E (E(^<^.J^^<^.0' 

+ E E (^"Ji ) ^ (^aj2 ) ^ (^aj3 ) ^ I 

<C(r2 + l)^^-l/2. 

At last, by the Cauchy-Schwarz inequality, (2.16), (4.3) and (4.54), we have 
IT23I <n-3EV2|E|(X-^X-^)_(ti)rtEV2|^ E(X-C/)„,(t2) 'I 



<n-/^EV4E(E(X„^^.)f }e^/^{E(E^.^.^«. 
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This completes the proof of the lemma. □ 

4.4. Central limit theorem for linear eigenvalue statistics of sample covari- 
ance matrices in the general case. Here we prove the CLT for the linear 
eigenvalue statistics of the sample covariance matrix not assuming that the 
fourth cumulant of its entries is zero [see (4.42)]. We use the scheme of the 
proof of Theorem 3.6 based on general differentiation formula (3.6) and the 
"a priory" bound (3.70) for the variance of statistics. Here is an analog of 
the bound for sample covariance matrices. 

Theorem 4.4. Let M be the sample covariance matrix (4-l)^(4-3) sat- 
isfying (4-40) CLnd (4-41)> ^6 corresponding truncated matrix (4-45), and 

Unr = Tr expjztM"^}. 

Then for any r > 0, 

(4.63) Var{n„,(t)} < a(/i4)(l + 1*1^)' 
and 

(4.64) VRr{Afnrm < CAf^i) (/ (1 + dt^ , 

where Ct{h4,) depends only on r and 

We omit the proof of Theorem 4.4, because it repeats with natural mod- 
ifications the proof of Theorem 3.5 for the Wigner case, and is again based 
on the use of the interpolation matrix (4.46) and known bound (4.29) for 
the Wishart matrix. 

Theorem 4.5. Let M he the sample covariance matrix (4-l)-(4-3) satis- 
fying (4-40) and (4-41) J and (/9:M^M he the test function satisfying (3.40). 
Then the centered linear eigenvalue statistic J\f° [<^] of M converges in dis- 
trihution, as m,n—> oo, m/n^ c^ [l,oo), to the Gaussian random variable 
of zero mean and variance 

(4.05) vsci.i = vw,.,.[.i + (/;; ,,)\ 

where Vwishiv'] is given hy (4.28), K4 = fi^ — 3a^ is the fourth cumulant of 
entries of X. 

Proof. Using the notation of the proof of Theorem 3.6, we note first 
that, according to Theorem 4.4 analogs of estimates (3.94) and (3.95), yield- 
ing the uniform boundedness of and dYnr/dx remain valid in this case. 
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To estimate dYnr/dt^ we write [cf. (3.96)] 

^ a,k 

with 

(4.66) $«fc = (y^C/^),fc(t)e°,(x), |E{Z)ifc^>,fe}|<CKt,x), / < 5, 
and by using (3.6) and (4.59)-(4.61), we obtain an analog of (3.98): 

at n , 

where, in view of (3.7), (4.61), (4.62) and (4.66), the error term is bounded 
by C2{t,x) in the hmit (4.7). The term Ti was calculated while deriving 
(4.33): 

Ti = ia^CnYnrixjt) + ia^tE{n~"'^n'„,-(t)e°^(x)} 

+ ia^ [ B{n~^Uj^^{t-ti)unT{ti)elr{^)}dti 
Jo 

— 2c?x J ti(^(ti)E{n~^u^^(t + ti)e„T-(x)} dti. 

We also have, by (2.16) and (4.3), 

B{\n-\Ut)\'}=n^'n\TTM^U^t)\'} 

(4.67) ^ 

<n-3E|^(^X;,.X;,) |<C 

I- j V a /J 

and, by integrating by parts, 

/ E{n"^'u'„^(t - ti)'u„^(ti)e°^(x)} dti 
Jo 

= f'-E{n~^u'^^{t-ti)}Ynr{x,ti)dti 
Jo 

+ /* E{<,(t - ti)n-i<,(ti)e°,(x)} dt 
Jo 



1) 



where the r.h.s. is uniformly bounded in view of (4.63) and (4.67). Hence, Ti 
is uniformly bounded for any r > 0, and so does d/dtYnr- This and analogs of 
(3.94) and (3.95) imply the existence of a subsequence {Ymr} that converges 
uniformly to a continuous Y^. 
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Now an argument similar to that leading to (3.104)-(3.106) yields an 
analog of (3.104): 

Ynr{x, t) = Tl, + rj^^„ + £^r,n{t^ x) + o(l), 

(4.68) 

?i,m— >oo,m/n— >c, 

where the terms on the r.h.s. are given by the r.h.s. of (3.105)-(3.107) with 
instead of w'^{l + 5jk), ^ak of (4.66), and for any r > the reminder term o(l) 
vanishes in the limit (4.7) uniformly on any compact set of {t > 0,x G M}. 

The term T\ was in fact calculated in the proof of Theorem 4.2 and is 
equal to Yn{x^t) of (4.34) and (4.35) with the Wishart matrix M replaced 
by the truncated sample covariance matrix M"^. Using (4.63) to estimate 
the reminder term r„ of (4.35), and noting that by an analog of (3.58) 
v^j- vmp ill the limit (4.7), we get an analog of (3.108) in the same limit: 

T:2^^ia\c-1) f'Yr{x,ti)dti 
Jo 



+ 



2ia^ [ VMp{t-ti)Yr{x,ti)dti + xZr{x)A{t) 
Jo 



with A{t) defined in (4.37). 

Consider now the term T^^ „ of (4.68), given by (3.106) with $ofe of 
(4.66). It follows from (4.50)-(4.53) and an argument similar to that of 
the proof of Lemma 4.1 that the contribution to „ due to any term of 

n~'^J2ak-^ak^'^i<^^ Containing at least one element {Y'^U"^)ak, vanishes as 
m,n ^ CO, m/n — > c. Thus, we are left with the terms, containing only diag- 
onal elements of and Y'^U'^Y'^'^ . These terms arise from enrD^f^{Y'^U'^)ak 
and 3L'afc(>^^t/^)afc^afce°^, and by (4.50)~(4.53), their contributions to rj^^„ 
are [cf. (3.110) and (3.111)] 

(4.69) 



a,k 



* [U^,, + i{Y^U^Y^' U * C/^J)(ii)e°,(x)} dh 
and 



a.k 



^E|e..(x)^ iU^,+iiY^U^Y^^U * U^k)iti)dti 

(4.70) 

X I h^{t2){Ulu + i{Y^U^Y^^)^^ * Ulk){t2) dt2 

Thus, the entries of C/^ and y^f/'^y^^ are present here in the form [cf. 
(3.112) and (3.113)] 

(4.71) KpO = E{VnT{h,t2)Wp^n{t3,ti)e°n^{x)}, p=0,l, 
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and 

(4.72) Kp = E{VnT{tl,t2)Wp^n{t3,U)enT{x)}, P=0, 1, 

where VnT{ii-,t2) is defined analogously to (3.114) and satisfies an analog of 
(3.115), and 

(4.73) Wp,n{hM)=n^^Y.(YUY-^)^o.{h){YUY-^)l^{U) 

a 

satisfies 




by (4.3) and (4.54). Since the expectations of VnT{ti,t2) and Wp^nitsjti) 
are uniformly bounded, and by Lemma 4.2 below their variances vanish in 
subsequent limit (4.43), then, applying the Schwarz inequality and (3.116), 
we conclude that 

Kpo = o(l), 

(4.74) 

Kp = ZnT{x)VnT{tl,t2)'Ei{Wp^n{t3,U)} + 0{l), P = 0, 1 

[cf. (3.118) and (3.119)], where the error terms vanish in the limit (4.7) 
uniformly in {t,x) G M?. Using the interpolation argument similar to that 
in the proof of Lemma 3.1 with the GOE matrix replaced by the Wishart 
matrix, we get an analog of (3.120): 

(4.75) lim VnT(tl,t2) =VMp{ti)vMp{t2)- 

m,n~*oo,m/n-~*c 

To find the limit of E{i(;o,n(i3i ^4)}, we note that 

where Vnrit) converges to VMp{t) as m,n ^ 00, m/n — > c, and that by (4.67) 
and a similar argument, the sequences {v'^ri^)} and {t'nr(0} are uniformly 
bounded, so that we have 

lim B{v'^^{t)}=i-Wj,jp{t) 

m,n—fOC,m/n—tc 

uniformly in t, varying in a finite interval. Furthermore, it can be shown by 
an argument, used not once before and based on (2.14), (3.6) and relation 
(4.80) below, that the functions 

lim E{{Y-U^{Y-f)^a{t)} 

and 

hm eUti-^ ^(y^c/-(y-)^),„(t)l = {ic)-\'Mp{t) 
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satisfy the integral equation 

h{t) = a^VMp{t) + ia^ / h{t — ti)vMp{ii) dh. 



This and Proposition 2.1(v) imply that the functions coincide, and we ob- 
tain, in view of (4.80), 

lim E{wi^n{t3,ti)} = -C~^v'Mp{t3)v'jyjp{ti). 

m,n^oo,m/n—*c 

We conclude from the above that the contribution of (4.69) to Ti^^^m vanishes 
as n; ^ c« uniformly in t and x, varying in any compact set of > 0, x S R}, 
while in (4.70) we can replace Ukk by vmp and {Y'^U'^Y'^'^)aa by {ic)~^v'j^p. 
We obtain 

(4.76) lim T:^^^ = -c-'K4xZr{x)C[v] f' A^,{ti)dti, 

m,n—iOO,m/n~^c JQ 

where 

(4.77) A^^{t)=cvMp{t)+ I VMp{t-ti)v'Mp{ti)dti, 

Jo 

(4.78) C[^]=iJt(p{t)A^,{t)dt 

or, in view of Proposition 2.1, (4.20), (4.21) and (4.38) [cf. (2.81)], 

A^M = ^£^ e'f^'^Aa^c-ifi-an.fdij. 
Plugging the last expression in (4.78) and integrating by parts, we get 

J a. yj^a^c - (/i - flm)"' 

This, (4.68) and (4.76) lead to the integral equation for Yr[x,t) [cf. (3.127) 
and (3.128)]: 

Yr{x,t)-ic?{c-\) \ Yr{x,ti)dti-2ia^ [ VMp{t - ti)Yr{x,ti) dti 
Jo Jo 

= -xZr{x) [Ait) + K^C-^C{ip\ A^^ih) dtl^ + £3r{t, x), 

where satisfies (3.127). 

Now, to finish the proof, we have to follow the part of the proof of Theorem 
3.6 after (3.127) to obtain (4.65). □ 
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Lemma 4.2. We have, under the conditions of Theorem 4-5 in the limit 

(IV, 

(4.79) Var{C/,-,(t)} = o(l), 

(4.80) Var{(y-C/^y^^),„(t)} = o(l). 

Proof. The proof of (4.79) repeats with natural modifications the one 
of an analogous assertion for the Wigner matrix (see Lemma 3.1). It is based 
on the interpolation procedure and following from the Poincare inequality 
(2.21) vahdity of (4.79) for Wishart matrices. 

To prove (4.80), we consider 

Vnr{ti,t2) = B{iY-U^{Y-f)Uti){Y^U^Y-f)Ut2)}, 

putting in an appropriate moment t2 = —ti to get \'aLr{{Y'^U"^ {Y'^)'^)aa{ti)} ■ 
We have, by (3.6) and (4.56), 

VnAti,t2)=n-y^Y.^{X^k{YW^)ak{ti){YW^Y^n^{t2)} 

k 

= J2 nDak{{Y^Unak{ti){Y^UAY-fr^^{t2))} + ei.. 



k 



where |eiT-| < ^iC2{t,x) by an argument similar to that used in (4.57) and, 
by (4.52) and (4.53), the sum on the r.h.s. is 

i rVr,r{tl-s)Vnr{s,t2)ds + B{Vnr{tl){Y^UAY^f ZM} 

Jo 

+ i B{v°At,- s){Y-U-Y-^)Us){Y^UAY-nAt2)}ds 
Jo 

+ itin-'Vnr{ti,t2) + 2n-'E{{Y^U^Y^^)aa{ti + t2)} 

+ 2in~^ [ Vnritl + s,t2 - s)ds. 
Jo 

It follows from (4.3) and (4.54) that T4r(ti,t2) and B{{Y^U^{Y^f)l^{t)} 
are uniformly bounded. This, the Schwarz inequality and (4.63) imply that 

E{vnAh){Y^UAYn^)Ut2)} = 0{n-') 

and 

|EK,(ti - s){Y^UAY^fU{s){Y^UAYn^ra^{t2)}\ 

2- 



<B(^\vUti-s)\(j2{Y^j, 
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j 

as m, n — > oo, m/n — > c. Besides, a bit tedious but routine calculations, simi- 
lar to those in the proof of Lemma 4.1, yield the boundedness of derivatives of 
ynT{ti,t2) for any r > 0. Thus, there exists a subsequence {mi,ni) such that 
the limit V{ti,t2) = lim.m,,ni^oo V^«i,T(ii) ^2) exists and satisfies the equation 

fti_ 

V{ti,t2)=i / VMp{ti - s)V{s,t2)ds. 
Jo 

Now Proposition 2.1 implies that V{ti,t2) = 0. This completes the proof of 
the lemma. □ 
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